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The  purpose  of  this  thesis  is  to  present  a  systematic  study  of  the 
type  of  noise  spectra  in  solids  or  solid  state  devices,  which  can  arise 
from  "transport  noise,"  which  stands  for  noise  due  to  combination  of 
causes,  such  as  diffusion  of  carriers,  generation  and  recombination  (g-r) 
processes  of  carriers  in  the  bulk  and  at  the  surface.  The  techniques  to 
be  described  apply  equally  well  to  the  problem  of  heat  diffusion,  with 
heat  transfer  or  reflection  at  the  boundaries  of  the  system;  we  have, 
however,  mainly  the  carrier  noise  problem  in  mind. 

After  an  introduction  we  describe  the  various  mathematical  techniques, 
applicable  to  transport  noise.  The  methods  studied  are  (1)  the  eigen- 
functions  expansion  method,   (2)  the  Green's  function  method.  Each  of  these 
methods  can  be  applied  to  two  physically  distinct  procedures,   (a)  the 
sources  or  Schottky  method,   (b)  the  correlation  method. 
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In  chapter  3  we  discuss  the  physical  noise  sources, viz.  the  diffusion 
source,  g-r  source  and  surface  source,  as  well  as  Clarke  and  Voss'  P- 
source . 

In  the  next  few  chapters  we  discuss:  transport  noise  in  symmetrical 
embedded  bodies,  in  nonsymmetrical  embedded  bodies,  in  symmetrical  non- 
embedded   bodies  and  in  nonsymmetrical  ncrsn4>edded  bodies.  The  concept 
"embedded"  body  refers  to  a  system  which  is  part  of  a  very  large  system, 
with  which  it  exchanges  particles  or  heat,  such  as  a  volume  between 
probes  on  a  large  semiconductor.  The  concept  "nonsmbedded  body"  refers 
to  a  body  inside  which  there  is  transport  noise,  with  deterministic 
boundary  conditions  on  the  surface  of  the  body. 

In  chapter  4  we  first  review  the  older  work  on  the  circular  patch 
and  the  cylinder  (MacFarlane,  Burgess,  Richardson,  van  Vliet,  Chenette) . 
Next  we  solve  for  the  noise  of  an  embedded  sphere,  caused  by  diffusion 
and  voltime  g-r  processes.  For  infinite  carrier  lifetime,  the  results  are 

found  to  be  in  agreement  with  those  of  Fassett.  For  finite  lifetime  an  ex- 

-3/2 
tension  of  his  results  is  obtained.  The  high  frequency  co     law  is 

confirmed. 

In  chapter  5  we  discuss  the  results  for  a  rectangular  embedded  bar, 

disc,  or  cube.  Because  of  the  lesser  symmetry,  the  spectrum  cannot 

be  obtained  in  closed  form,  but  only  in  a  triple  integral  form. 

The  result  has  been  numerically  evaluated  by  computer  techniques.  Not 

all  the  ranges  predicted  by  Clarke  and  Voss  are  visible,  however. 

For  the  P-source  the  noise  has  also  been  evaluated.  Contrary  to  Clarke 
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and  Voss'  claim,  no  1/f  range  is  visible.  Thus,  their  heat  diffusion 
noise  model  cannot  account  for  the  observed  1/f  noise  in  metal  films. 

In  chapter  6  syiranetrical  imembedded  bodies  are  discussed.  Parti- 
cular attention  is  given  to  the  noise  in  a  sphere,  with  boundary  condi- 
tions due  to  surface  recombination.  A  closed  form  is  obtained  via 
an  expansion  of  the  Green's  function  in  spherical  polynomials. 

In  chapter  7  the  nonsymmetrical  bar  is  discussed.  The  results 
given  previously  by  Lax  and  Mengert  and  others  were  evaluated  by 
computer. 

In  chapter  8,  finally,  we  discuss  noise  due  to  the  surface  g-r 
source.  This  noise  can  be  dominant  due  to  the  surface  modulation 
factor;  it  has  not  been  considered  in  all  older  theories.  We  discuss 
three  exan^les:  (a)  a  semiconductor  single  crystal  sample  with  surface 
g-r  noise;  Cb)  the  noise  in  the  channel  of  a  MOSFET,  caused  by  diffusion 
into  the  adjacent  bulk  and  by  g-r  (surface)  processes  near  the  oxide 
layer;  Cc)  the  same  model  but  with  special  care  for  the  alteration 
of  tLe  diffusion  constant  and  mobility  in  the  channel  due  to  surface 
scattering.  For  models  (h)   and  Cc)  it  is  shown  that  a  realistic  large 
1/f  range  (10  hz  up  to  10   hz)  can  occur.  Thus  this  model  can  account 
for  much  of  the  high  frequency  1/f  noise  in  MOS  devices. 

The  other  exanples  and  models  considered  in  this  thesis  firmly 
indicate  that  volume  g-r  and  diffusion  processes  give  rather  smooth 

spectra,  but  without  a  1/f  range.  The  high  frequency  asymptote  is 

-3/2  2 

always  co     (embedded  case)  or  co   (nonembedded  case). 


IX 


CHAPTER  1 
INTRODUCTION 


1.  1   Object 

The  purpose  of  this  study  is  to  make  a  systematical  investiga- 
tion of  diffusion  and  generation-recombination  noise  spectra.  Such 
studies  may  aid  in  understanding  the  physical  mechanisms  underlying 
'l/f  noise.  The  noise  spectra  obtainable  from  transport  processes 
are  analyzed  and  the  possibility  of  contributions  to  'l/f  noise  is 
studied.  It  will  be  shown  that  most  transport  processes  caused  by 
volume  sources  do  not  give  rise  to  l/f  noise. 

In  addition,  a  detailed  investigation  has  been  undertaken 
of  surface  noise.  This  is  noise  caused  by  surface  sources  which 
give  rise  to  stochastic  boundary  conditions  at  the  surface.  Surface 
noise  contributions  are  studied  in  detail,  since  the  older  theories 
due  to  Hyde^,  Champlin^,  Lax-Mengert3,and  van  Vliet**  only  considered 
the  volume  noise  contribution.  This  analysis  gives  for  the  first 
time  a  transport  noise  model  which  results  in  l/f  noise;  however, 
the  frequency  range  for  which  it  occurs  is  often   too  high  to  ac- 
count for  the  experimental  data,  unless  very  long  lifetimes  of  the 
carriers  can  be  shown  to  exist. 

1 .2   Flicker  Noise 
The  phenomenon  of  l/f  noise  or  'flicker  noise'  started  with 
the  study  of  J. B. Johnson^  on  noise  in  vacuum  tubes  in  1925.  He 
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noted,  that  besides  shot  noise,  which  has  a  white  frequency  spectrum, 
at  low  frequencies  another  noise  component  appears,  which  varies 
roughly  inversely  proportional  to  the  frequency.   The  first  interpre- 
tation of  this  noise  was  given  by  W.Schottky^  in  1926.   Schottky 
attributed  the  effect  to  fluctuations  in  the  work  function  of  the 
cathode  material,  caused  by  the  random  arrival  and  leaving  of  foreign 
atoms  at  the  cathode  surface.  The  spectrum  computed  by  Schottky 

was  not  exactly  of  the  form  S(f)  =  A/f;  rather,  his  effect  results 

2  2 
m  a  Lorentzian  spectrum,  S(f)  =  B/(l+to  t  ). 

Since  Schottky 's  mechanism  failed  to  interpret  Johnson's  1/f 

data,  as  well  as  the  flicker  noise  later  on  observed  in  a  large  number 

of  solids  and  devices,  various  authors  have  proposed  a  great  variety 

of  detailed  studies  to  understand  the  physical  mechanisms  behind  1/f 

noise.   In  the  subsequent  chapters  we  shall  summarize  some  of  these 

efforts  and  extensively  study  the  contribution  of  transport  and 

surface  noise  to  the  origin  of  flicker  noise. 

1.3   Transport  Noise 
Transport  noise  is  caused  by  quantities  subject  to  stochastic 
transport  in  the  sample.   The  most  notable  examples  involve  particle 
diffusion  and  heat  diffusion.   Both  are  described  by  a  linear  partial 
differential  equation.  The  noise  can  be  obtained  in  two  ways.   In 
the  first  procedure,  the  'sources'  or  'Langevin'  procedure,  one  adds 


noise  source  terms  to  the  phenomenological  transport  equation,  which 
then  becomes  of  the  form 

(|^+  A)aCi:,t)  =  ?C£,t)  (1.3.1) 

2 
where  A  =  -DV   and  a  =  AT  or  An  for  heat  or  particle  diffusion 

respectively,  and  ?  is  the  volume  noise  source.  The  spectrum  of  the 
latter  must  be  known. 

In  the  correlation  procedure,  on  the  other  hand,  one  uses 
the  ensemble- averaged  stochastic  equation,  which  by  Onsager's^  prin- 
ciple, is  the  phenomenological  equation.  Thus  we  must  then  solve 

C|^+  A)  <aCr,t)|aCr,0)>  =  0  (1.3.2) 

where  ^|^  denotes  the  ensemble  average  of  a(r,t),  conditional 
to  a  given  value  a(r,0)  at  time  t  =  0.   In  this  case  the  covari- 
ance  function 

riX'X')  =  <AaCr,t)Aa(;£',t)>  (1.3.3) 

must  be  known. 

The  connection  between  both  procedures  is  given  by  A -theorem 
due  to  van  Vliet^,  which  connects  the  source  spectra  to  the  cova- 
riance  function 

CA^.-^A^)rCr,r')  =  ^^Cr,r').  (1.3.4) 

To  solve  either   (1.3.1)   or   (1.3.2)   one  can  use  two  mathematical 
formalisms,   viz, the  eigenfunction  expansion  method  and  the  Green's 


function  description.  Thus,  in  total,  there  are  four  methods  to 
obtain  transport  noise.  The  equivalence  of  all  these  methods  has 
been  shown  by  van  Vliet  and  Fassett'» *".  The  details  of  these 
methods  are  described  in  the  next  chapter. 

1.4   Surface  Noise 
This  is  noise  caused  by  stochastic  boundary  conditions  at 
the  surface.  An  example  occurs  when  a  semiconductor  sample  is 
subject  to  surface  generation-recombination,  with  the  surface  re- 
combination velocity  S  being  unequal  to  zero  or  infinity. 

In  the  experiments  on  modulation  of  conductance  by  surface 
charges,  it  was  found  that  only  10  percent  of  the  induced  charges 
were  effective  in  changing  the  conductance.  J.Bardeen^^  proposed 
that  the  ineffective  portion  of  the  induced  charges  is  lodged  in 
states  localized  at  the  surface.  Since  the  work  of  McWhorter^^ 
and  others,  it  is  held  that  the  traps  in  the  oxide  modulate  the 
generation-recombination  rate  and  this  gives  rise  to  flicker 
noise  in  the  device.   Kooge^^  explained  the  excess  noise  by  an 
empirical  model,  though  in  his  case  the  noise  is  a  volume  rather 
than  a  surface  effect. 

Recently  van  Vliet,  van  der  Ziel  and  Schmidt  I'*  considered 
heat  conduction  in  a  thin  film  supported  by  a  substrate  and  found 
that  surface  sources  can  give  rise  to  1/f  noise.  A  sample  with  a 
small  inversion    layer  near  the  surface  which  is  coupled  to  larger 


bulk  volume  is  considered  in  this  thesis.  The  model  is  mathematically 
similar  to  that  considered  by  van  Vliet  et  al.  >   and,  as  for  the  thin  film 
case,  the  surface  noise  gives  rise  to  1/f  type  of  noise  behavior   with 
some  theoretical  restrictions. 


1.5   Embedded  and  Non-embedded  Bodies 


From  the  theory  it  emerges  that  one  must  distinguish  between 

models  involving  embedded  and  non-embedded  domains.   In  the  former 

case,  one  supposes  that  the  volume  of  interest  for  the  noise,  V  ,  is 

part  of  a  larger  domain  V.   Fluctuations  in  the  particle  number  or  heat 

content  of  V  arise  due  to  sinks  and  sources  in  V  (g-r  contributions) 
^  s 

and  due  to  transport  across  the  boundary  into  the  rest  of  V  (diffusion- 
drift  contributions) .  An  example  is  afforded  by  the  case  of  field 
emission  from  a  metal  tip.  Here  ions  within  a,  say  circular,  patch 
of  the  active  area  tip  can  modulate  the  electron  emission;  these  ions 
therefore  affect  the  emission,  when  they  are  within  this  circular 
patch  V  ,  which  is  a  part  of  the  larger  cathode  area  V  over  which 
ions  diffuse. 

The  non-embedded  case  applies  when  we  monitor  the  particle  or 
heat  content  of  a  domain  V  which  does  not  communicate  with  the  ex- 
terior but  on  the  boundary  of  which  we  have  'deterministic  boundary 
conditions',  usually  of  the  Dirichlet,  Neumann  or  mixed  type.  An 
example  is  when  we  measure  noise  due  to  carrier  fluctuations  in  a 
semiconductor  sample  V^,  subject  to  surface  generation-recombination 
(mixed  boundary  conditions)  on  the  free  surface. 


This  dissertation  is  divided  into  three  parts.  The  first 
part.  A,  contains  the  mathematical  and  physical  basis;  in  particular 
in  chapters  2  and  3  we  discuss  the  various  mathematical  methods  for 
transport  noise  and  the  physical  noise  sources,  respectively.  The 
second  part,  B,  gives  the  noise  due  to  volume  sources.   In  chapter  4 
we  discuss  the  noise  for  symmetrical  embedded  bodies.  Chapter  5 
is  devoted  to  noise  from  nonsymmetrical  embedded  bodies.  Chapters 
6  and  7  consider  the  noise  from  symmetrical  and  nonsymmetrical  non- 
embedded  bodies,  respectively.  Finally,  in  the  third  and  last  part, 
C,  we  consider  the  noise  due  to  stochastic  boundary  conditions. 
In  chapter  8  we  study  in  detail  the  noise  caused  by  surface  generation- 
recombination  processes  in  bulk  samples  as  well  as  in  MOSFETs. 


PART  A 
MATHEMATICAL  METHODS  AND  PHYSICAL  NOISE  SOURCES 


CHAPTER  2 
MATHEMATICAL  METHODS  FOR  TRANSPORT  NOISE 


2.1  Methods 


As  indicated  in  the  introduction,  there  are  basically  two 
ways  to  attack  the  problem  of  transport  noise.  The  first  class  of 
procedures  considers  (a)  transport  equation(s)  in  the  sense  of 
Langevin,  that  is,  one  adds  source  terms  ?(r, t)  which  themselves 
have  no  memory,  to  the  phenomenological  equation Cs)  for  the  fluc- 
tuating variableCs)  aC^,t);  the  variables  a(r^,t)  are  Markovian 
random  of  infinite  dimensionality  since  ^  is  a  continuous  parameter 
like  a  position  in  some  domain  V.  The  boundary  conditions  on  the 
surface  S  of  V  may  be  deterministic  or  stochastic;  in  the  latter 
case  surface  Langevin  sources  ^(r  ,t)  must  be  specified  (r  €  S) , 
The  "Langevin  procedure"  or  "sources  procedure"  is  still  the  most 
versatile  one.  The  stochastic  Langevin  equation  can  be  solved  with 
Green's  functions,  orthogonal  or  biorthogonal  expansions  or  Fourier 
methods . 

The  second  class  of  procedures  is  based  on  the  solution  of 
the  transport  equation  for  the  ensemble  averaged  transport  quanti- 
ties <a(r,t)  |a(r,0)>  =  ^^^'^^^afr  1'  ^"  which  the  fluctuations 
Aa(^,t=0)  are  specified.  The  solution  involves  the  "correlation 
procedure."  The  noise  is  obtained  by  the  Fourier- Lap lace  trans- 
formed Green's  function  of  the  transport  equation,  providing  the 
covariance  function,  TC^,^')  =  (Aa{;^,t) AaC;^' ,t)>  is  known; 


'  i-'ii': 


here  (  >  denotes  a  stationary  non conditional  average.  This  leads 
to  the  form  suggested  by  van  Vliet  and  Fassett^.  This  procedure 
is  mathematically  much  faster  than  the  Langevin  procedure;  however, 
it  has  the  disadvantage  that  stochastic  boundary  conditions  cannot 
be  included.  Nevertheless,  if  the  transport  equation  is  linear  and 
if  the  surface  sources  are  uncorrelated  with  the  volume  sources, 
the  former  can  be  separately  considered  with  the  sources  method,  due 
to  the  superposition  principle.  Also,  despite  the  fact  that  the 
covariance  function  in  the  presence  of  surface  sources  is  often  not 
unique,  it  has  been  shown  by  van  Vliet^°  that  the  spectra  resulting 
from  the  Green's  function  procedure  are  unique.  We  will  discuss  both 
classes  of  methods  and  their  prospects  for  obtaining  1/f  noise. 

2.2   Langevin  Equation  and  Source  Methods 
The  stochastic  transport  equation  is  supposed  to  be  of  the  form 

LaCr,t)  =  C|^+  ApaCr,t)  =  ?(x,t).  (2.2.1) 

where  L  is  a  linear  differential  or  integral  operator.  For  the  case 
of  diffusion: 

a  =  An  Cparticle  diffusion) 
A  =   -D7   ,  (2.2.2) 

a  =  AT  (heat  diffusion) 

In  the  diffusion  model  A  is  a  self  adjoint  operator  which  simplifies 
the  mathematics  considerably.   The  boundary  conditions  must  be 
specified  at  the  boundary  S  of  V.  The  noise  may  concern  the  entire 
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domain  V  or  subdomain  V  c  V.   If  V  is  infinite,  the  only  condition  is 
that  a(x,t)  goes  sufficiently  fast  to  zero.   If  on  the  surface  S  of  V 
we  have 

"'aC^,t)  =  0,  x,Q   ^   S,  C2.2.3) 

m  being  a  linear  surface  operator,  the  boundary  conditions  are  deter- 
ministic; if 

maC^,t)  =  K(z^,t),       Xo  ^  S,  (2.2.4) 

the  boundary  conditions  are  stochastic. 

2.2.1  The  Bulk  Noise 

Because  of  the  superposition  principle,  noise  associated  with 
^  and  ^  can  be  separately  considered.  Thus  first  setting  ^  =  0, 
the  boundary  conditions  are  deterministic  (2.2.3).  Various  methods 
will  be  given  to  solve  from  (2.2.1)- (2.2. 3) . 

2 . 2 . 1 (a)  Eigenfunction  expansion 

Let  4>^(r)   be  the  eigenfunctions  of  A  subject  to  b.c.  (2.2.3) 
that  is 

U^(Ji)   =   Vk^'  (2.2.5) 

1®^  ^k^-^  ^®  ^^®  eigenfunctions  of  the  adjoint  operator  A^,  then 

A\(r}  =  X*tj^(r).  (2.2.6) 
The  functions  4>   and  \|^  are  biorthogonal,  i.e. 

C«j,,V  =/  ^'^\(Z)n(T)    =6^^.  (2.2.7) 


E. 
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Expanding     a(x>^)   and  K(x>^)   i"  ^^^   functions   (L    and  denoting  the 
expansion  coefficients  by  a,  (t)   and  \W  ,  t:hen   (2.2.1)   gives 

2[^Ct)+\aj^(t)]*j^(^   =  2  ^it)4>^00-  C2.2.8) 

k  k 

Multiplying  with  i»(r)  and  integrating  over  V  and  with  (2.2.7)  one 
obtains 

C^+  \)a^(t)  =  ?^(t),  (2.2.9) 

which  is  a  standard  one  variable  Langevin  equation.   In  the  usual 

method  of  solution  truncated  Fourier  series  in  the  time  interval 

CO,T)  are  used.  Fourier  integrals  are  not  suitable  since  the 

Langevin  equation  (2.2.9)  does  not  apply  for  negative  T;  thus 

00      ico  t 
a^(t)  =  Z  a^^e  ",00^  =  2nn/T,  n=0,±l,...,     (2.2.10) 
n=-oo 

and  similarly  for  ^j,(t),  one  obtains 

<.   at  >    =  Jfknf^n^s 

^  Kn  ^n'S        ,-•     ^   ^  ^   •     ^*^ 

Now  the  essence  of  the  sources  method  is  that  the  white  spectra 
of  (?CG,t)CCE' ,t)),  denoted  as  S„C£,^'),  are  known.  Thus  according 
to  the  Wiener- Khintchine  theorem,  the  sources  correlation  function  is 

<?Cr,t)?*(r',f)>^  =  ls^Cr,^')5(t-f)  (2.2.12) 

CC=C*  is  real),  in  which  S  must  be  known.  Expanding  in  E.F.,  we 
have 


K- 
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Z  (^(t)C*(f  )>^<^^Ct)<i>*(r')  =  ^  (x,r)6(t-f ).       C2.2.13) 

3  3 
Multiplying  this  with  ^^j^,  (x)i/i  Cg')  >  integrating  over  d  rd  r' , 

and  using  biorthogonality,  it  is  found  that 


Representing  the  C's  by  truncated  Fourier  series  we  further  have, 

i(Go  t-oo  t')  ico  (t-f) 

nm 

(2.2.15) 

from  which  we  obtain 

<VL>s  -  2T^//d'rd\'S^(X.r)i:C?:)t,Cr').      (2.2.16) 
The  spectrum  is  given  by 


\a*^V  =   '^^'^<\n^ln>s-  ^2.2.17) 


Thus  from  (2.2.11)  and  (2.2.16) 

S.  .*(«)  =  r-  ^^  ^i  • — T7T //d"^rd^r'S^fr,r')>J/*rr)A|r.(r').(2.2.18) 

Finally  for  the  spectra  of  ^a(r,t)a(^' ,t))  one  obtains  for  the  noise 
due  to  the  bulk  source  ^: 

Sa"^''Cr,r',co)  =  Z  S   *(co)*^Cr)<^*(r'),  (2.2.19) 

K  V    K  V 

or  substituting  from  (2.2.18)  one  obtains 

Qbulk,   ,   ,   ,  4>^Or)4>l(r')        ^^  ^3  ,3  ,  ,     ,,*..,.. 
^   ^'^  '   ~~  U    (\-i^)  C^|-i")  -^^   ''l  ^2ScCri,r2)ti^Cri)i,(r2)  . 

(2.2.20) 
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To  show  the  connection  with  the  correlation  procedure  to  be 
discussed  in  section  2.3,  we  write 

(y ico)  (\*-iM)      =      XjJT^  (jy^  ^   Xjli^j    •  (2.2.21) 

For  the  covariance   from   (2.2.20)    and    (2.2.21) 

~~     2  L  "^VT"^^"  'l'  ^2\Ctl'^2^^k(^PnCr2). 

(2.2.22) 
Operating  on  the  LHS  with  A  +A  ,  we  note  with  (2.2.5) 

(2.2.23) 
Observing  the  closure  property 

2  ctj^(r)>^*(r^)  =  5(r-rp  (2.2.24a) 

K 

2  </>*(r')i^(r2)  =  6(r'-r2)  (2.2.24b) 

and  integrating  out  the  delta  functions  we  find 

(Y^^-^^^"^''^'^'^  =  |«gCr.r').  (2.2.25) 

This  is  the  A-theorem  due  to  van  Vliet^.  This  theorem  gives  the 
connection  between  the  sources  spectra  and  the  covariance  function. 

Also  it  warns  that  sources  spectra  cannot  be  postulated  ad  hoc,  if 

„bulk,   ,^  .   , 

r    (X'Z  J    IS  known  from  statistical  mechanics. 


^• 
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Substituting  (2.2.25)  in  equation  (2.2.20)  and  using  equation 
(2.2.21)  one  obtains 

.bulk,   ,   ,    ,  ,  \ix)n^n        1      1  ^ 

X  ffd^r^d\(A^^A     )T^'''^^U^,j^-)-^*,Or,^n^iX^).    (2.2.26) 

Introducing  Green's   theorem  for  the  operators  A  and  A    : 

(Af,g)-(f,A+g)    =  #C[f,g*].da,  (2.2.27) 

where  ^C-da     is  the  bilinear  concomitant  involving  a  surface  integral. 
Thus  in  (2.2.26) 

^^d^r^C[T^'''^^(T^,T^),^^(T^)].  (2.2.28) 

Doing  the  same  for  At^ ,     and  keeping  in  mind  the  closure  property, 
we  get 

Sa"^''Cr,r',co)   =  2  Z^^Cr)  ^^^i^ /d^r"rCr",r' )^|/*  (r") 

K         K. 

+  hcj  +  surface  terms,  (2,2.29) 
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where  hcj   stands  for  Hermitian  conjugate,   r  -»•  r' ,   ico  -»■  -Im, 
X  ->■  X*;  the  surface  term  involves  the  bilinear  concomitant.  The 
first  part  of  (2.2.29)  is  called  the  volume  contribution.   It  will 

be  denoted  by  S   .  Wlien  the  bilinear  concomitant  is  zero  -  as  is 

■^      a 

often  the  case  -  the  volume  contribution  equals  the  bulk  contribu- 
tion [2. 2.  29). 

Often  the  interesting  part  is  to  know  the  spectra  of  the  total 

value  or  spatial  average  of  variable  a(r,t)  in  V  .   Thus  let 

ACt)  =   /  aCr,t)d\.  (2.2.30) 

V 
s 

The  integration  yields  for  (2.2.20) 
Sources  form: 

.bulk, ,     ,  ,  ,  ,3 ,3 .  h^'^n^'^ 

S.        (m)   =  Z  /  /  d  rd  r'  -px .  .  ,,  ..  .  s 

^  k/v  V  (Xj^^ico)(X|-ico) 

s  s 

X  /  /  d^r^d\2S^(Xl>X2)ikC£l)^^^U2^  (2.2.31) 

3 
The   correlation   form,   obtainable   from   (2.2.29)    contains  /d  r'r(^",^'). 

Let  us  suppose,    as   is  often  the   case, 

rCt",^')    =    [<AA^>/V^]5C?:"-r').  (2.2.32) 

Then  we   find 


-liiM 
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Correlation  form: 

sy°Vco)      =     ^\,     ^  Re  /     /     d\d\'    !•     \      .     C2.2.33) 

^  ^s  V     V  k      \'^"' 

s     s 

In  general   X,    =    Mv  +  iv,      is   complex;      for  self  adjoint  case     v,    =  0 

and   (2.2.33)   reduces  to 

vol  Ubkh        \\^'^^^\^lO\^ 

Svol(co)      =     ^V  ^   -^     ^     ^  •  (2.2.34) 

a  V,22 

s       k        M, +00 

Finally,  note  that  the  eigenfunction  method  has  several 
drawbacks.   If  the  domain  V  is  finite,  the  results  occur  in  a  series  form 
which  is  not  easily  evaluated.   If  the  eigenfunctions  can  be  taken 
to  be  plane  waves,  then  we  end  up  with  the  integrals  over  Jk- space 
which  are  not  trivial.  We  therefore  prefer  the  Green's  function 
method,  whenever  the  Green's  function  can  be  evaluated  in  closed 
form,  as  is  the  case  for  suitable  geometries  (one  dimensional  bar, 
cylinder, sphere) .  The  limiting  spectral  ranges  can  then  be  obtained. 

2.2.1(b)  Green's  function  procedure 

The  Green's  function  procedure  is  much  more  direct.   The 
Langevin's  equation  (2.2.1)  is  Fourier  transformed  with 

icO  t       jrt-n 

aC^.t)  =  Z  a^C?:)  e  "  ,«^=-^,  n=0,±l,±2, .  .  .      (2.2.35a) 
n 


ioj  t 
?(X,t)  =  2  ?^Cr)  e  '^   .  (2.2.35b) 

n 


F" 
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This  gives 

Cico^+A^)a^(X)    =  ^n'^^-  (2.2.36) 

The  Green's  function     G(x,X,' >i"  )      is  defined  by 

(ico^+A^)GU,X'.%^    =   5C?:-X'),  (2.2.37) 

and  satisfies  homogeneous  boundary  conditions  corresponding  to  equa- 
tion (2.2.3).   Using  (2.2.37),  the  solution  of  (2.2.36)  is  found 
by  using  the  standard  procedure 

a^(r)  =  /  d^^G(r,ri,ia3„)^(rp.  (2.2.38) 

The  spectiumis  found  by  the  result   analogous   to    (2.2.17) 

Sa"^^;^.^'.")   =   lin>  2T<a^(x)a*(;r')> 
T-x» 

=  /  /  d^rjd^r2G(;^,^^,ico)G(r',r2,-ico)S^(r^,r2);  (2.2.39) 

for  S,         we  find 
A 

S^"^^(«)   =  /     /     d^rd^rV  /  d^rjd\2^(X,r^,ico)G(r',r2,-ia.)S^(rj,r2). 

^      ^  (2.2.40) 

These  are  the  sources  Green's  function  forms. 

It  is  well  known  that  Green's   functions   can  be  expanded  in  a 
biorthogonal  series 

*kC£)K(r') 
GCr,r',i&i)      =     2  -A: A .  (2.2.41) 

^*~    '  ,  \,  +100  ^  ^ 


f 
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When  (2.2.41)  is  substituted  in  (2.2.39),  the  result  (2.2.20)  is 
recovered. 

2.2.2   Surface  Noise 

Now  we  consider  the  transport  equation  (2.2.1)  with  ^  =  0 
and  while  we  have  stochastic  boundary  conditions  (2.2.4);  we  also 

assume  that  ?  and  K   are  uncorrelated.   (If  this  is  not  the  case  there 

cross 
will  be  a  contribution  S     .) 

3. 

2. 2, 2 (a)  Eigenfunction  expansion 

This  method  has  not  been  tried  due  to  the  complexity  of  a  series 
expansion  in  eigenfunctions;  however  a  response  method  without  Green's 
functions  can  be  derived. 

2.2.2(b)  Green's  function  procedure 

Again,  we  make  a  Fourier  analysis  of  the  (homogeneous)  Langevin 
equation  and  of  the  stochastic  boundary  condition  (2.2.4) 

(ico  +A  )a  fr)  =  0,  (2.2.42) 

^  n  _;^'^  n^~"^ 

"'^^^=^n^,^-  ^'-'-''^ 

The  solution  of  these  equations  leads  to 

a    (x)   =  fda  ti(T,T   ,iw  )K   fe  )  (2.2.44) 

n^^      0    ^^l~o   n  n^^o 

2 
(da  =  d  r  ),  where  H  is  the  surface  Green's  function.   For  the 
0       0    ' 

spectra  we  now  obtain  due  to  the  surface  sources 
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S^^^'^CCXS")   =  f^^%  da'H(x,X^,ico)H(x',;c;.-i")S^U:o';^).  (2.2.45) 

while  for     S^"'^       one  obtains 
A 

Sr^'^Cco)    =  /     /     d\d\'#  ^  da  da  H(x,r   ^coDHCt'.Xo'-i'^^S    (r   ,x:). 
^  V     V  S  S 

s     ^  C2.2.46) 

In  particular,    let     S^  =  ZCCo)^'^^^  (^-^)  •     '"'®" 

Sf'^^c.)  =  #  iV^^I  ^  d\H(x,r^,ico)|^         C2.2.47) 

O  V 

s 
These  are  the  sources  Green's  function  forms.  The  covariance 

function  associated  with  the  surface  noise  is 

r-^'(X,X')  =   /  ~Re  S-^^Ct,^',.).  C2.2.48) 

surf 
It  can  also  be  shown  that  T         (x>/G')  satisfies  the  homogeneous 


A-theorem 


(A^+A^,)  r'^^'^CC.r')  =  0   (r,r'  €  V).  (2.2.49) 


'Ny    'X/ 


htillc      sinrf 
Total  covariance     T  =  T         +r         ,   therefore, satisfies  the  A- theorem 

which  is 

(A^+A^,)r(x,x')    =  ^^(X,X')    =  S(x,r'),  (2.2.50) 

where  H  is  the  correlation  strength  of  the  volume  sources. 

2.3   The  Phenomeno logical  Equation  and  Correlation  Methods 
The  phenomenological  equation  corresponding  to  (2.2.1)  is 

^<^Cr,t)>^^^,3  =  (|^  ^V<^Cr.t)>^^^03  =  0,       (2.3.1) 
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where  ()  ,  ^,.  means  a  conditional  ensemble  average  for  systems 
with  fixed  fluctuations  at  time  t=0.  The  boundary  conditions  cor- 
responding to  (2.2.3)  or  (2.2.4)  is 


m 


<^f^'^»a(;C,0)  =  '■  ^'-'-'^ 


2,3.1  The  Volume  Noise 

Since  the  surface  sources  average  to  zero,  it  is  unlikely 
that  the  noise  found  by  correlation  method  includes  surface  noise. 
Moreover,  the  correlation  method  does  not  give  the  complete  bulk 
noise;  in  fact,  it  only  gives  what  is  identified  as  the  volume  noise 
in  section  2.1.  The  reason  for  this  is  that  the  covarismce  func- 
tion r,  which  figures  centrally  in  this  method,  only  represents  the 
solution  of  the  inhomogeneous  A-theorem.   The  surface  source  and 
other  boundary  spectra  give   covariances  which  satisfy  the  homo- 
geneous A-theorem. 

2.3.1(a)   Eigenfunction  expansion 

The  stochastic  variable  is  once  more  expanded  in  eigenfunc- 
tions  (i>,{T),   which  satisfy   equation  (2.3.2),  this  yields  for 
for  equation  (2.3.1) 

<\f^^>{a^(0)}*\<^kf^^>{aj^(0)}  =  0'  f2.3.3) 

which  gives 

<\^^)>{a^(0)}=<\f°)^^  f2.3.3a) 

and  furthermore 

\,it)    .   <a^(t)a|(0)>^  =  «\f^)>{a^(0)}^|tO^>  s  =  ^  'V\CO)a|(0)>^ . 

(2.3.3b) 
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Finally  for  the  solution  one  obtains 


Ca       =  2  /dt[e-^"\a^(t)a*CO)>^  +  e^"^<a^(0)a|Ct)>] 
k   V  0 

=   2<a,a*V     (^-^+  ^Att).  C2.3.4) 


■<V!>s    ^T^^ip^^ 


We  have  used  here  the  transposition  property  for  the  correlation 
function  ^,p(,t)   =  "^o,  (-t).  With  an  expression  similar  to  (2.2.19) 
we  find 

SI°'U.Z'^-)   -   ?J  *kte)*^*(X')<VPsf^  ^  Xjk^-  ^2.3.5) 

Next  we  consider  the  expansion  for  the  covariance  function 
rC£i.X2)  -   <AaC£^,t)Aa(X2,t)> 

=  1   <\^Vsh^l^n^2^  -  (2.3.6) 

Multiplying  with  l^j^CXi)  '*!'';  (Xo)  ^nd  integrating  over  all  r^.    and  ^_, 
one  obtains 

<a^,a*,>   =  ffd\^d\TOi^:ir^n^,(z^n^,(T^).  (2.3.7) 

Substituting  ttis  into  equation  (2.3.5)  we  obtain 

K  V  V  V 

Carrying  out  the  required  integration,  the  result  is  obtained  as 

^l°^(Z>Z''^)    =2  Z  /dV'c#,^(x)ij^(r")(^;-i^)r(r",r')  +  hcj  ,      (2.3.9) 
k  V  k 


^' 
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which  is  the  same  result  as  obtained  previously  via  the  sources  method 

3 
in  equation  (2.2.29).   For  the  result  of  the  noise  of  A(t)  =  /  d  r  a(^,t), 

V 

using  this  relation  one  obtains  (2.2.33)  or  (2.2.34).        ^ 


2.3.1(b)      Green's   function  procedure 

For  this  procedure  we  need  the  solution  of  equation  (2.3.1); 
for  this  purpose  we  Fourier- Laplace  transform  equation  (2.3.1),  i.e. 
let 

00 

£(^,i")     =     /  e"^"^<aC£,t)>^^^Q^   dt.  (2.3.10) 

Then  equation  (2.3.1)  yields 

Cico+A^)£(X,iM)  =  a(x,0).  (2.3.11) 

The  Green's  function  is  given  by 

Cico+A^)GC5,X',ico)  =  5(x-X'),  (2.3.12) 

with  the  boundary  condition  mG(^  ,^',ia5)  =  0  corresponding  to 
equation  (2.3.2).  By  the  standard  procedure  we  find  the  solution  as 

aCt.i")  =  /  d\"G(x>X".i'^)a(x",0).  (2.3.13) 

V 

Therefore,  by  inversion  of  the  Fourier- Laplace  transform 

<^^'^^>aCr,0)  =  ■  k^'^'^^''^   /d^r"G(x,^",i^)a(x",0)   (2.3.14) 

where  co  is  a  complex  frequency  and  C  is  contour  encircling  all  poles. 
For  the  correlation  function  we  have  (t  >  0) 
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=   1   ^  d(5e^"^  /  A"G(x,X."4")r^"^^C?:',X").      (2.3.15) 
^   C       V 

For  negative  time  (-t)  we  use  the  transposition  property 
<^CG>X'>-'t)  =  't>(x'  >Z>^)  •  For  ■the  spectrum  the  Wiener-Khintchine 
theorem  then  yields, 

rjVol,        ,      .  1    ^j.j  J-   i(co-(Jo)t    ^j3  ,,„,       ,,    .-^„bulk,    ,      ,,^       ,     . 

^a     l^^'I;  '")   =   -  ^  ^^^^  ^^  ^^  r"G(x,X  ,i")r         (r',r")   +  hcj  . 

0     C  V 

(2.3.16) 
Interchanging  the  time  and  contour  integral  we  have 

/  d^/  dt  e^^'^-"^^G(x,^".iw)   =  -^  d^[^^i-r-+tT5(^-co)]G(x,X"4") 
CO  C  il-w-wj 

=  -2TrGC5,^",ico+0)(by  Cauchy's  theorem)  (2.3.17) 

So  now  substituting  the  integral  in  equation  (2.3.16)  we  get  the 
result 

S^°^(X»X'.")  =  2  /d\"G(x,X",i'^)r^"^''(X'«X")  +  hcj.  (2.3.18) 

V 

This  is  van-Vliet-Fassett  form.  The  result  is  in  full  accord 
with  equation  (2.3.9)  if  we  use  the  series  form  of  the  Green's 
function,  equation  (2.2.41). 

2.3.2  The  Boundary  Noise 

The  spectrum  found  by  the  above  procedure  does  not  include 
surface  noise.  However,  it  also  does  not  comprise  all  the  bulk  noise. 
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which  is  to  some  extent  surprising.  The  reason  is  that  the  transport 
equation  is  not  Markovian  in  the  usual  sense.  The  ;fact  that  the  di- 
mensionality is  infinite  means  that  the  usual  Markov  property 

W2(aC;C,t),a'(x',f))   =  P(aC;c.t)  |a' (x,' ,t'))Wj(a' (r' ,f ))  ,    (2.3.19) 

which  is  the  basis  for  the  evaluation  of  the  correlation  function  in 
(2.3.15)  must  be  doubted.  In  other  words,  the  process  cannot  be 
Markovian  unless  boundary  behaviour  is  also  specified. 

The  total  noise  can  be  broken  up  in  two  ways .  The  first  way 

is 

gtotal  ^  gbulk  ^  gsurface  ^^gcross^ _  (2.3.20) 

This  breaking  up  is  according  to  the  considerations  of  section  2, 
the  terms  representing  effects  due  to  the  bulk  sources  C(r,t),  the 
surface  sources  ^(^,t)  and  due  to  their  possible  correlation.  In 
the  second  way  we  evaluate  volume  noise  according  to  section  (2.3.1) 
based  on  (2.3.15).  The  remaining  part  to  be  evaluated  in  this  section 
is  called  s^°""dary   ^^^ 

gtotal  ^  gvolume  ^  gboundary^^gCross^_  (2.3.21) 

„       ,^    .       „boiindary  ^     ^      -^u  ^u 

To  obtain     S  -^     we  start  with  the  response  expression 

(2.2.39).      In  the  RHS  we  substitute   for  S„  the  A-theorem,   and  obtain 

S^"^^(X,r'.»)   =   2  /  /d^rjd\2^Cr,rj,i«)GCr',r     -ico)(A^+Ar)rCri,r^). 
V  V  'vl  ~2 

(2.3.22) 
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We  indicate  by  f(jr^  ,^)  here  any   (interior  domain)   particular  solution 
of  the  A  theorem,  which  may  differ  from  physical  covariance  function 
by  a  solution  of  the  homogeneous  A  theorem.     The  RHS  of  C2.3.22)   is 
rewritten  as 

RHS  =  2  /  d\     G(j^;e,i.i")    /  dV  G(^' ,^,-ico)  (-icof  A^  )f(r,,r  ) 

V  -^V  i.2 

+  2  /  d^^  GixJ  ,J^.-io:>)   f  d\  GCr,r^,ico)  (;  cofA^  )fCr^,r2) 

V  V  ~1 

=  2  /  d\     G(x,Xi,i«){/  d\  fCr     r  )L^     G(t' .r     -ico) 

V  Vi2 

4^   d2rQC[fC;;^,^Q),GC5>,^Q,-ia))]} 

+2/  dVcC^'^     -ia)){/  dVfCr  ,r Jl];*  GCr,r     i^) 

V  V~l 

+^   d^rQClfC5Q,^2)'GCr,rQ,ioo)]},  (2.3.23) 

where  we  have  used  Green's  theorem  in  the  form 

CL  f,g)-(f,  I'^g)  =  (Af,g)-(f/^g)  =#dj  C[f,g*].  (2.3.24) 

According  to  the  definitional  equation  for  the  Green's  function 

L^  GCE',^2'-i'^)  "  ^^  ^  C52''5''^"^   =  5Cr2-r'),  (2.3.25a) 

L  *  GC5,^^,ico)   =  L  *  G*Cr  ,r,ico)   =  6  Cr  -r)  .  (2.3.25b) 

"^l  ~1 

Carrying  out  the  delta  function  integrations  we  are  left  with 

S^"^''C5>5',co)      =     2  /  d\"  GCr,r",ico)fCr",r')    +  hcj 

V  " 

+2/dV  GC£,r",ico)  #  d^r^  C[f  Cr",rQ)  ,G(r' ,r^,-ico)]   +  hcj  .    (2.3.26) 
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The  first  two  terms  comprise  S   .  The  last  two  terms, 

a. 

together  with  S^"^^  comprise  5^°""^^^^^^.  The  Green's  function 
form  for  S^"^   was  given  in  (2.2.45).  We  then  have 

S^""'^^.^'^  "3   =  2/  d\"  GC?:,r",ic.)  #  d\   C[fCr",Xo),G(x',Xo'-i"3] 
^  V  s 

+  hcj  +  /  /  d\  d\'   f  f   d^r  d^r-  HCt,Xo,ico)H(x',Xo'-i")-  (2.3.27) 

V  V  S  S 

s  s 

We  note  that  all  terms  in  (2.3.27)  are  quadratic  in  the  Green's 

m.    ij-.-    /7f--^  nVol   J  „bound 
function.  Thus  the  division  of  S  into  S    and  S       separates 

the  spectra  in  expressions  linear  and  quadratic  in  G.  The  merit  of 
this  separation  is  that  in  many  instances  all  terms  of  S      ^  are 
zero  for  homogeneous  systems  and  delta  covariance  function.  For  in- 
homogeneous  systems  such  as  transitors,  the  boundary  terms  are  of  great 
importance,  howeverjwhen   S  °""     =0,  the  evaluation  of  the  noise 
by  S^°   as  given  in  section  (2.3.1(b))  is  by  far  the  fastest  pro- 
cedure . 

In  equation  (2.3.26)  we  used  a  particular  solution  T   .  Another 

'^  vol 

choice  of  this  solution,  i.e.,  F' ,  leads  to  another  division  of  S 

and  s  °""  ''y^.  However,  the  total  result  (2.3.26)  is  unique.  Usually 

one  chooses  T   such  that  the  spectrum  S        vanishes,  if  possible. 

If  on  physical  grounds  it  can  be  argued  that  there  is  no  boundary 

A 

spectrum  contribution,  then  this  choice  of  T   is  clearly  the  physical 
covariance  function. 


27 


So  far  we  have  discussed  the  mathematical  methods  to  obtain 
transport  noise  spectra.   In  the  next  chapter  we  shall  discuss  and 
briefly  explain  the  different  physical  noise  sources,  their  spectra 
and  importance  in  calculating  overall  noise  spectra  of  the  physical 
system. 


CHAPTER  3 
PHYSICAL  NOISE  SOURCES 


3.1   Introduction 
We  have  seen  in  the  previous  chapter  that  a  calculation  of 
noise  spectra  with  the  Langevin  or  the  sources  method  requires  that 
the  spectra  of  the  Langevin  noise  source  must  be  known .  There  are 
various  physical  noise  sources  with  known  spectra.  Furthermore,  we  have 
also  seen  that  there  are  two  different  types  of  noise,  viz.  first 
volume  or  bulk  noise  and  second   surface  or  boundary  noise  stemming 
from  volume  noise  sources  and  surface  noise  sources,  respectively. 
In  this  chapter  we  shall  try  to  explain  both  types  of  noise  sources 
as  we  encounter  them  in  physical  systems. 

3.2   Volume  Noise  Sources-^ ^ 

3.2.1  Particle  Diffusion  Noise  Source 

To  explain  the  diffusion  noise  source  let  us  assume  that  the 
fluctuations  in  number  of  particles  (electrons  or  holes)  are  governed 
by  the  diffusion  equation,  say  for  electrons 

~r   CAn)  ♦  DV^CAn)  =  0  C3.2.1) 

ot 

where  An  is  the  fluctuation  in  density  of  electrons. 
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The  Langevin  noise  source  representation  of  this  diffusion  equation 
in  the  form  of  the  standard  diffusion  noise  source  ?(x>t)   is  given  by 

1^  (An)  +  DV^CAn)  =  ?(x,t) .  (3.2.2) 

The  source  C  follows  from  the  stochastic  current  as  follows.  For 
(AjJ)  ,. -_  we  have  the  randomness 

<^^ji^diff  "  -qDV(An)  +  qr)(x,t)  .  (3.2.3) 

On  a  microscopic  scale,  "nCtjt)  represents  the  thermal  scattering  or 
Brownian  motion  velocity  fluctuations  of  the  carriers.  This  r\  source 
has  a  spectrum^ ^ 

S   =  4D<TiC5)>6C5-^')I  (3.2.4) 

where  I  is  the  unit  tensor. 

We  also  write  <n(^)  =  n^(jc)  .     The  conservation  theorem  for 
charge  reads 

^  -  -  V-(AJ)  ....  =  0.  (3.2.4a) 

9t   q  ~   ~  diff 

Substituting  (3.2.3)  we  obtain 

^+  DV^(An)  =  V.TiC£,t).  (3.2.5) 

Comparing  (3.2.5)  with  C3.2.2),  we  note 

^(X't)    =   v-Tic?:,t). 
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Hence  for  the  spectrum 

S^U.X')  =  2^^,    ■■   S^Ct.X')  (3.2.6a) 

=  4D2^-S;^,{nQ(x)5(x-X')>.  (3.2.6b) 

This  diffusion  noise  source  is  the  proper  noise  source  of  sta- 
tistical mechanics.  As  we  can  see  by  the  A-theorem,  due  to  van  Vliet^, 
which  connects  Langevin's  noise  source  spectra  to  the  covariance 
function, 

^V^X'^^^"^'^   =  |s^U:.X')  (3.2.7) 

the  standard  diffusion  noise  source  gives  for  the  covariance  function 
a  delta  function: 

rix.X':)  =  nn(r)5(r-r«).  (3.2.8) 

All  the  previous  calculations  of  diffusion  noise  by  McFarlane^*, 
Burgess     ,   Richardson  ^,   van  Vliet  and  Fassett^,   Lax  and  Mengert^.are 
based  on  the  assun^tion     that   the     spatial  covariance  function  is  a  lielta 
function  given  as     constant  x  5(^-^').     This  result  can  also  be  proven 
from  first  principles,  see  ref.9,  page  329. 

3.2.2     Heat  Diffusion  Noise  Source 

The  heat  diffusion  noise  source  can  be  explained  by  considering 
the  relation  between  the  heat  conductivity  a  and  the  heat  current   K. 
which  is 

K     =      -aVT  +  rjIX.t)  (3.2.9) 
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where  "H  is  a  Langevin  source.  Heat  conservation  requires  that 

cd|f  +  V-K,  =  0  (3.2.10) 

where  c  is  specific  heat  per  gram,  and  d  is  the  density.  We  can 
find  the  stochastic  heat  diffusion  equation  by  substituting  (3.2.8)  into 
(3.2.9)  and  by  using  T  =  <T)+AT, 


3  AT 


aV 


2(AT)   =  "^^^S^   =  5CE,t)  (3.2.11) 


3t    ^  ''      ca: 

where   a  =  a/cd  is  the  thermal  diffusion  constant.   The  correlation 
strength  or  white  noise  spectrum,  of  the  source  ri  has  been  found  by 
van  Vliet?  using  the  Boltzmann  transport  equation,  to  be 

§nfe'^   =  4Ik  T(^  ^o(T)5(r-Tn,  (3.2.12) 

where   T(x)   is  the  time  or  ensemble  averaged  temperature  at  position 
X,   and  I  is  the  unit  tensor.   Clearly 

S^Ct,^')      =      C4k/cV)    V^.V^,{  T(x)    ^a(r)5Cr-r')}  (3.2.13) 

which  is   the  spectrum  of  the  standard  heat  diffusion  noise  source.^'* 

3.2.3     Voss   and  Clarke's  P-Source^^ 

Voss   and  Clarke  introduced  another  noise  source,  which  has  a 
delta  function  type  correlation  function  and  a  white  spectrum  given  by 

S^v-ci    ^'^'^      =     Po^Cr-r')  (3.2.14) 

where  P     is   a  constant. 
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As  we  have  seen  previously  noise  sources  cannot  be  chosen  ad  hoc 

due  to  limitations   imposed  by  the  A-theorem  and  statistical  mechanics. 

Using  the  A-theorem  the  three   dimensional  spatial   covariance   function 

corresponding  to  the  P  noise  source  is   found  to  be 

p2 
r(x,X')     =     j^    ^^y|    .  C3.2.15) 

Such  a  long-range  correlation  function  is  highly  unphysical. 


3.2.4  Particle  Volume  Gene rat ion -Re comb in at ion  Source 

Usually  noise  due  to  generation  and  recombination  of  carriers 
results  in  fluctuations  in  the  resistance  of  the  specimen  under 
consideration.  The  noise  can  be  described  by  the  fluctuating  number 
N  of  the  carriers,  say  for  electrons;  then  the  equation  governing 
the  fluctuations  in  N  is 

^  =  g(N)  -  r(N)  +  Ag(t)  -  Ar(t)  (3.2.16) 

where     gCN)   and  r(N)   are  generation  and  recombination  rates  and  Ag(t) 
and  ArCt)    are  the  randomness   in  these  rates. 

Substituting     N  =  N  +AN,   where     N_     is  the  equilibrium  number 
of  carriers   and  neglecting  higher-order  terms   in  AN     yields  the 
equilibrium  condition     g(N/v)    =  r(N^)      and  the   linearized  Langevin 
equation,       given  by 

^    .     -f  .48(t)    -   4r(t)    ,   -  f  .  C^_,(t)  (3.2.17) 

where 

T  MN       dN^|N      • 
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The  Langevin  source  C  (t)   can  be  obtained  by  assigning  shot 
noise  to  the  random  transitions  Ag  and  Ar: 

S,g  =   2g(NQ)  .  2gQ  (3.2.18) 

hr     -     2r(No)  =  2rQ  =  2go  (3.2.19) 

S,   ,   =0    (no  correlation). 

Ag ,  Ar 

Hence 

Sr    =  S,  +  S.   -  2S.   .   =  4g„.  (3.2.20) 

'^gT  Ag    Ar     AgjAr    "0 

Making  a  Fourier  analysis  of  (3.2.17)  one  finds  for  the  spectral  connec- 
tion 2 

S  -r^  4g  T 

^  (I+CO'^'C  )       1+W  T 

Integrating  over  all  frequencies,  we  find 

which  is  Burgess'2°>2^  g-r-theorem. 

In  transport  systems  the  spatial  dependence  must  be  added. 
Thus,  for  the  density  nC5,t)   subject  to  diffusion  and  generation 
recombination  processes  we  have 


dAn 
8t 


^.  DV2  An(x,t)  -  ^2i£^    =  ^diff^'^^  ^V^''^ 


The  g-r  noise  source  now  becomes 

S.   Ct,X')  =  4g5(^-r').  (3.2.22b) 

gr 
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3.3  Surface  Noise  Sources 


3.3.1     Particle  Surface  Generation-Recombination  Source 

Just  like  generation-recombination  processes  in  the  bulk  give 

rise  to  a  bulk  Langevin  noise  source,  so  generation-recombination 

processes  at  ths  surface  can  be  represented  by  a  surface  noise  source; 

this  latter  source  causes  surface  transport  noise,  at  the  surface 

of  the  sample. 

For  example,    for  the  case  of  non-zero  surface  recombination 

velocity,  the  stochastic  boundary  condition  can  be  written  as 

D^V^CAn)  +  Sto  =  ^Uo^t)  (3.3.1) 

where  D     is  thediffusion  constant   for  electrons         V„     is  the  surface 
n  '"'O 

concentration  gradient  and  ^(^,t)    is   the  surface  stochastic  Langevin 
noise  source.      For  the  threedimensional   case  this  surface  noise 
source  has   a  spectrum 

S^     =     4MaSn(rQ)5f2^Cr-r^)  (3.3.2) 

and  for  the  one  dimensional   case 

S        =     4MaSn(XQ)/A 

where  M  is  the  modulation  factor  >  1,  S  is  the  surface  recombination 
velocity,  A  is  the  area  of  surface,   a  =  <AN  )/N   and  5'-^-'c^-^) 
is  the  two  dimensional  delta  function. 
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3.3.2  Heat  Transfer  Source 

The  heat  transfer  source  is  a  stochastic  source  due  to  ab- 
sorption and  emission  of  blackbody  radiation  at  the  surface  of  a 
body.  For  example,  this  occurs  in  the  case  of  a  thin  film  supported 
by  a  substrate  with  the  front  face  of  the  body  being  blackened, 
so  that  it  sees  the  blackbody  radiation  field.  The  boundary  condition 
in  this  case  becomes 

3  AT 
c^CO)   ^  -   gATCO)      =     Ut).  (3.3.3) 

'x=0 
Here  ^(t)   is  the  stochastic  heat  transfer  source.     The  spectrum  of 
this  source  is   given  by  the  generalized  Nyquist's   theorem 

\     =     4kT2g     =     16agkT^  (3.3.4) 

3 
where     g  =  4a  T       is  the  heat  transfer  conductance,  T     is  the  ten^era- 

ture  of  the  environment  and  the  mean  temperature  of  the  body,     and 

5  4         3  2 
Og  =  2n  k  /15h  c       is  the  Stefan-Boltzmann  constant. 


PART  B 
NOISE  DUE  TO  VOLUME  NOISE  SOURCES 


CHAPTER  4 
SYMMETRICAL  EMBEDDED  BODIES 


4.1  Introduction 
In  this  chapter  we  shall  study  the  noise  spectra  resulting  from 
symmetrical  geometries  which  are  part  of  an  infinite  expanse  i.e.  the 
"embedded"  case.  We  have  already  seen  in  a  previous  chapter  that  there 
are  two  volume  noise  sources,  namely  the  standard  diffusion  noise 
source  and  Voss  and  Clarke's  'P'  source.  Further  we  have  also  seen 
that  the  'P'  source  is  highly  unphysical.  In  this  chapter  we  con- 
sider only  noise  stemming  from  the  diffusion  noise  source. 

4.2  Survey  of  Older  Work 

In  this  section  we  briefly  summarize  the  results  of  previous 
work  for  obtaining  noise  spectra  from  symmetrical  embedded  bodies 
like  the  linear  bar,  the  cylinder  and  the  sphere,  prior  to  the  dis- 
covery of  the  Green's  function  method. 

MacFarlane   in  1950,  studied  the  problem  of  the  power  density 
spectra  of  noise  current  resulting  from  contact  noise  in  semiconductors 
or  from  emission  noise  in  field  emitter  diodes.   In  these  theories 
low  frequency  noise  is  attributed  to  the  random  movement  of  absorbed 
ions  on  the  contact  surface  between  two  granules  of  a  carbon  resistor 
or  on  the  emission  tip  of  a  field  emitter  tube.  The  current  through 
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the  tube  or  the  emission  of  electrons  is  assumed  to  take  place  only  at 
localized  patches  on  the  surface  and  the  absorbed  ions  are  assumed  to  give 
rise  to  a  dipole  layer,  which  modulates  the  work  function  in  such 
a  way  that  AW  is  linearly  related  to  the  concentration  of  ions  in 
that  patch.  This  concentration  fluctuation  is  due  to  diffusion  of 
ions  over  the  surface  of  which  the  conducting  patch  is  a  small  part. 

MacFarlane  found,  by  first  calculating  the  autocorrelation 
function  and  then  using  the  Wiener-Khintchine  theorem,  that  the 
spectral  power  density  of  the  noise  current  is  given  by 

RH  -  /{l-expC-x^/'^)[lQ(x"^/Vli(x'^/'*)]}cos(px)  dx       (4.2.1) 

2 
where     p  =  cox  ,     t     =    (2r)   /D,   r  is   the  radius  of  the  patch,   and 

D  is  the  diffusion  constant  of  ions,   and  I's  are  the  modified  Bessel 

functions.      In  the   limiting  cases  it  was  shown  that,  for  a  circular 

patch,    (using  probability  theory  results      .) 

R(M)      oc     0.2345  p"^^^,     p  «   1  (4.2.2a) 


and 


R(co)      cc     0.1466  p"^^^,     p  »  1  (4.2.2b) 


while  for  a  long  thin  strip 


and 


RCco)      cc      (2p)"-^/^,  p  «   1  (4.2.3a) 


R(w)      cc      (2p)'^/^,  p  »  1.  C4.2.3b) 
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For  the  long  thin  strip  this  shows  the  occurrence  of  the  3/2  power  law 
at  high  frequencies.  This  universal  law  was  violated  by  the  result 
C4.2.2b);  Burgess^ ^showed  MacFarlane's  result  to  be  in  error;  when 
properly  integrated,  a  3/2  power  law  at  high  frequencies  also  result 
for  the  circular  patch. 

4.3  Previous  Work  Using  the  Green's  Function  Method 

4.3.1  Linear  Bar 

Van  Vliet-Fassett^  considered  similar  problems,  as  discussed 
above,  using  the  Green's  function  procedure.   For  the  one-dimensional 
diffusion  equation 

^     -     oi^.  (4.3.1) 

ax 

where   Ap  is  a  fluctuation  variable  and  D  is  a  diffusion  constant, 
the  Green's  function  is  found  by  standard  procedure  to  be 

GCx,iw,x')      =     -^^  expCiUglx-x' I)  C4.3.2) 

where     u^  =  Vco/D  e  . 

The  spectrum  for  the  fluctuations  AP,  where  P  is  the  total 
number  of  particles  in  a  segment  L  is  found  to  be 

2 

SapC"^)  =  -^^O-  Cl-e'®(cose+sine))  (4.3.3) 

^         DG 


^ 


40 


where  9  =  L(-^)   ,  in  agreement  with  Burgess'  solution ^^     of  the 


.16 


integral  given  by  MacFarlane  .  The  low  frequency  spectrum  goes  as 

-1/2  -3/2 

CO     and  the  high  frequency  asymptote  is  again  (Xi         .     The  latter 

behavior  is  characteristic  for  all  diffusion  spectra. 


4.3.2  Cylinder  or  Circular  Patch 

Van  Vliet  and  Chenette^^  studied  the  noise  spectra  resulting 
from  diffusion  processes  in  cylindrical  geometries.  The  spectral 
densities  are  derived  from  the  relevant  Green's  functions.  Also 
extensions  to  MacFarlane' s^^  spectrum  for  the  circular  patch  were 
obtained. 

The  noise  spectrum  for  the  fluctuating  variable  X,  where  X 
is  the  total  number  of  particles  in  the  cylindrical  volume  V  ,  in 
terms  of  the  relevant  Green's  function,  is  given  by 

S^Cco)  =  ^^^f     f     GC?:,ico,r')dv  dv'  (4.3.4) 

s   D  D 
s  s 

where  V     is  the  volume  of  domain  D     (cylinder  or  circle),  which  is 

part  of  an  expanse  V  »  V   . 

The  Green's  function  was  found  to  be  the  Hankel   function  of 

the  first  kind  of  order  zero;   the  spectrum  is  given  by 

„   f  -.  8n  CvarX) a^  _         1        ,1         i„(l)  r     ■. ,   .     -.-, 

S^Co.)      =  D  ^^  7— T  f-  27  *  2^1   ^Cxa)J    (xa)], 

(xa) 

(4.3.5) 

where  H  and  J  are  Hankel  and  Bessel  functions,  a  is  the 

radius  of  the  cylinder,  D  is  the  diffusion  constant  and  x  =  iV^^^^  . 
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It  was  found  that  for  low  frequencies  the  spectrum  has  loga- 

-3/2 
rithmic  behavior  and  at  high  frequencies  the  co    behavior  takes 

over. 


4.3.3  Sphere 

The  problem  of  three  dimensional  diffusion  has  been  solved  by 

van  Vliet  and  Fassett^.  Using  the  method  of  reciprocity,  the 

2 
spectriim  was  found  to  be  (for  fluctuating  variable  AP,  with  AP  =  Pq) 

I^^Pq'^O   2    -r  2 
G.p(co)  =  p {y  -2+e  [y  CcosY+siny)  +4ycosy  +2(cosY-sinY)] }, 

Dy 

^  C4.3.6) 


where 


Pq     -     average  number  of  particles  in  the  sphere 
R^     -     radius  of  sphere 
D       -     diffusion  constant 


and 


Y     =     2RqVCw/2D) 


At   low  and  high  frequencies  the  spectrum  was  found  to  have  the  limiting 

values : 

32TrP  R^ 


PQC2D)^/24TrR2 


Again,   at  high  frequencies  the  spectrum  follows  the    'universal' 
3/2  power  law. 
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4.4  Present  Work:  Carrier  Fluctuations  in  a  Sphere 
Due  to  Diffusion  and  Generation-Recombination 

In  this  section  we  investigate  the  noise  spectra  resulting  from 
fluctuations  in  the  carriers  in  a  sphere,  due  to  diffusion  across  the 
boundaries  and  due  to  sinks  or  sources  inside  the  sphere.  This  problem 
has  not  been  solved  hitherto. 

We  consider  a  sphere  of  radius  'a'  (volume  V  )  which  is  part 
of  a  very  large  expanse.  The  fluctuating  variable  is  x,  the  concentra- 
tion of  particles  in  the  spherical  volume  under  consideration;  the  total 
number  of  particles  in  the  sphere  is  X.  The  boundary  conditions  on 
the  surface  D  of  the  sphere  are  'fictitious'  i.e.  x     and  its 
derivatives  are  continuous  across  the  boundary. 

We  consider  the  equation 

|i  +  ax  -  DV  x  =  0  C4.4.1) 

ot 

where  a       =  t     is  the  volume   lifetime  and  D  is  a  diffusion  constant. 
The  Green's  function     g(i^,t  |i^' ,t ')   satisfies  the  equation 

|fCr>t|r',t')   +  agCr,t|r',t')    -  DV^gCr^t |r' ,t' )   =  5(t-t')6Cr-r') . 

(4.4.2) 

Assuming  t'  =  0   as  initial  time,  with  g(r,0 |r' ,0)=0,  and  taking  the 

Laplace  transform  of  equation  (4.3.2)  we  obtain* 

gGCr,s,r')  +  aGCr,s,r')  -  DV^GCr,s,r')  =  6(x-r').    (4.4.3) 

This  can  be  written  as 

'n^U:,s,x')  -  ^U,s,x')  +  v2G(x,s,r')  =  -  isCr-r')  (4.4.4) 
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or 


where 


(x^+v2)GC?:,s,x')   =  -^(Z-Zn  (4-4.5) 


=  i,^  ,  i^^  .  C4.4.6) 


Since  the  boundary  conditions  on  U  are  continuity  of  x  and  its 

derivative  and  the  domain  D  is  assumed  to  be  very  large,  the  ultimate 

boundary  condition  on  the  boundary  of  D  is  that  G  should  behave 

properly  at  infinity,  i.e.  go  sufficiently  fast  to  zero. 

Adopting  a  spherical  coordinate  system  we  have 

2 

r  r  sin  6  Bcp    r  sine 

Following  Morse  and  Feshbach  ,  page  808,  one  can  show  that  for  the 
infinite  domain  the  Green's  functions  in  spherical  coordinates 
only  depends  on  |r-r,'  |  =  R.  So,  evaluating  V  in  spherical  coordinates 
with  R  =  0  as  origin,  we  find 

^   -  ;!  aR  '-'^  8R^ 

K 
oK 

Also  we  can  write 

GCr,s,r')  =  G^CR).  (4.4.9) 

So,   the  Green's   function  equation   (4.4.5)  becomes 

^'^x^'^)        2   ^^xW  2,    ,,,  1   5(R)  ,,,,,, 
2~  *   R  ~9R~  "   ""  ^x^^^      =      "D          2"-  (4.4.10) 
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Following  Morse  and  Feshbach^^,  page  809,  we  need  a  solution  for  the  above 
equation  which  goes  as 

^  for  R  ->  0.  (4.4.11) 

Therefore  we  require  the  solution  to  satisfy  the  following  conditions 

i)      G  (R)  -> -^  as  R  -^  0.  (4.4.12a) 

In  addition,  the  noise  spectrum  for  high  frequencies  co  ->•  ■»   should  go 

to  zero,  or 

ii)        lim   G^(R)  =  0.  (4.4.12b) 

X  -*■  -00+ i" 

The  only  solution  for  equation  (4.4.10),  satisfying  the  above 

two  conditions  is 

ixR 
G^(R)  =  1^  .  C4.4.13) 

The  proof  of  the  solution  is  as  follows. 

^         .   ±ixR    ±ixR     , 

hx""  =  "tIrD- *  T- Ir  4  •  W.4.14) 

^2  o    T   ±ixR  ^   1     ±ixR  ^2   , 

9  /-  /-^>^     2„  ^  2ixe     6   1    e     9    1  ^a    .    -.r-^ 

-2G^(R)  =  -XG±— 5 __.__^__  C4.4.15) 

^   ^  T   ±ixR   -  ±ixR  ^   1 

RiR^x^'^^  -  *— ^2—   -DR—  9R  4^  •  (.4.4 A6) 

4nR  D 

The  second  term  of  (4.4.15)  can  be  shown  to  not  contain  a  delta  function; 
it  cancels  the  first  term  of  (4.4.16).  The  third  term  of  (4.4.15)  and 
the  second  term  of  (4.4.16)  are  combined  to  yield 

±ixR    2   7  .   1      ^±ixR  „  , 

h^^-^'m^^    -     ^T-'  4^  •  C4.4.17) 

oK 

For  the  latter  we  have 
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e^  V^  JL    =     -  V-  5C€-^')    =   -  ^^  •  C4.4.18) 

D  4TrR  D  ^-^^  4^yj^2 

To  prove   (4.4,18)  we  must  show 

R^V^  i  =  0     if     R  ^  0.  (4.4.19) 

K 

and 

/      V^  y^^    dh     =     -4Tr.  (4.4.20) 

small  sphere    '    ' 

The  first  statement  follows  from  substitution  of  (4.4.8).  For  the  latter 

statement  we  use  Gansz'  theorem: 

,i^'-F^^'^  =  ^^-^"fW  dS.  (4.4.21) 

Now    grad(l/|r-r' I)  =  -1/R   and  dS  =  4nR  ,  which  proves  (4.4.20). 

We  notice  that  the  solutions 

±ixR  '■ 

^^^5  =  ^4» 

satisfy  the  equation 

CV^+x2)G^(R)   =  -isCt-X')- 

Using  the  condition  of  equations (4.4.12a)  and  (4.4.12b),  we  find 

that  the  condition  (4.4.12a)  is  satisfied  by  both  solutions  but  using 

equation  (4.4.12b)  we  find 

±ixR      ±R(-i«'-") 
l^"*  .  ^4Wrd  =      4nRD  ^^.4.22) 

and  the  condition  is  satisfied  only  by  the  solution 
+ixR 

^x^^^  =  hw-  C4.4.23) 

which  is  the  required  solution  for  the  Green's  function. 
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For  the  fluctuating  variable  X,  the  noise  spectrum  is  given 
by 

S„C")  =  ^^^^^/  /  G(r,iw,r')  dV  dV  (4.4.24) 

s   D  D    "    " 
s  s 

The  required  double  integration  sought  for  the  noise  spectrum  is 

I  =  /  /  d\  d^r'  G  (R).  (4.4.25) 

V  V  "^ 

s  s 

Further  R  can  be  written   as 

R  =  v'r^+r'^-2rr'cose  (4.4.26) 

where  9  is  the  angle  between  r  and  r' .  Using  polar  coordinates  with 

r'  as  polar  axis  and  integrating  (4.4.25)  first  over  r  we  need  the 

integral  (denoting  by  I.) 

fl       2 
ixVr  +r'  -2rr'cos9 

(4.4.27) 

4TrDVr^+r'^-2rr'cose 

Further  we  write  this  integral  as 


.     [2      .2 
ixvr  +r 


-                                   ^^v^    ^^'    -2rr'cos9 
I^     =     -/r  dr  /d(p/d(cos9) zZZ=ZIIZi:r—  •  (4.4.28) 


)^A^•2 


4nDVr  +r'    -2rr'cos9 
Realizing  that  integration  over  cp  contributes  2Tr  and  further  that 


'  2       2 
,    ixVr  +r'    -2rr*cos9 
-ixrr'e 

I  2        2 
d                  ixvr  +r'    -2rr'cos9 

d(cos9)   ^ 

> 

/  2        2 
Vr  +r'    -2rr'cos9 

(4.4.29) 

we 

find  that    (4.4.28)    reduces  to 

T         -              ^            r    rdr    rr^'^^^*^'^ 

e  ^^^'-^  ].                       (4.4.30) 

h     -     2xD  r'i  {  ''^''   LC 
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As  to  the  double  integration  (4.4.25),  we  can  write 


a  a  a    r'  a    a 

/  /  f(r,r')dr  dr'  =  /  dr'  /  dr  f(r,r')  +  /  dr'  /  dr  f(r,r') 

0  0  0     0  Or' 

and  it  can  be  shown  that  the  second  contribution  equals  the  first  one,  i.e. 


a  a 


/  /  dr  dr'  f(r,r')  =  2  /  dr'  /  dr  f(r,r') 
0  0  0     0 


(4.4.31) 


because  of  the  symmetry  of  the  integrand  f(r,r').  We  can  thus  restrict 
ourselves  to  r  S  r'   and  omit  the  absolute  sign  in  (4.4.30). 

Realizing  the  fact  that  integration  over  the  solid  angle  dS2  con- 
tributes 4Tr,  the  required  double  integration  yields 

T    8T1  ^   ,  ,  , 
I  =  '^  ^■.   J    r'dr' 


2xDi 


ix(r+r')    ix(r+r')     ix(r'-r)    ix(r-r') 
re         e         re      "^e      ^ 

: + +  : _  


IX 


IX 


or 


(4.4.32) 


4Tr 

xDi 


/  r'dr'  [(f^.  i.)e2i^^'  -.  (f^-^)] 


X 


'IX    2' 

X 


(4.4.33) 


Using  the  following  relations  due  to  integration  by  parts 


/r'dr'  e 
and 


2ixr' 


r'e 


2ixr'    2ixr' 


2ix 


4x' 


-    T   ,      ,2  2ixr'    ,  2ixr' 
,2,  ,  2ixr'     r'  e       r'e        1 
'  dr'e      =  — —. +  = — 


/r'  dr'e 


2eixr' 


2ix 


2x 


4ix 


(4.4.34) 


(4.4.35) 
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the  required  integration  is  finally  found  to  be 


4n 
xDi 


2  2ixa     2ixa    2ixa     3     2 
a  e      ae      e       a     a 

+  :: —  +  ; + 


-2       .3     _4     3ix^2.4 
2x       IX      2x       '^   2x    2x 


(4.4.35) 


Using  for  the  noise  spectrum  equation  (4.4.24)  and  representing  the 

2 

AD       p  _ 

fluctuating  variable  as  p_,  with  —rr-  =  y     =   Pf^,  the  number  of  holes 

s    s 
per  linit  volume,  we  get 


S^p(co)=Re^l::P^  [encp2-4cp-.4).(\  -   cp2.4)  |  (4.4.36) 


5  p    ^         3 
Of 

where    cp  =  2ixa, 


Equation  (4.4.36)  can  be  rewritten  as 

64Trp  a  „    «   3 

S^p(co)=Re 1—  I4e'''-4cpe%^e%  +  \  "  4] .  (4.4.37) 

Dcp 

For  the  required  noise  spectrum  we  need  to  evaluate  the  real 
part  of  expression  (4.4.37).  We  shall  discuss  two  different  cases. 


4.4.1  Pure  Diffusion  (a=0) 

If  no  sinks  or  sources  are  present  then  a=0.  In  this  case 
the  required  expression  for  the  noise  spectrum,  taking  the  real  part 
of  equation  (4.4.37),  becomes 

16Tra  Pq  2      -y     2 

S.p(co)     =     g {-2+Y  +e      [r   (sinY+cosY)+4YCOSY+2(cosY-sinY)]  }, 

Dy 


(4.4.38) 


where       y  =     2aV^ 
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The  result    (4.4.38)    is  the  same  as  previously  obtained  by 

van  Vliet  and  Fassett^  for  the  case  a=0.     The   low  and  high  frequency 

asymptotes  are  found  to  be   (expanding  sin  and  cos  terms) 

32np-a 
^Apf")       I.-15D-  (4.4.39a) 

"^^  PoC2D)l/24na2 

^Apf^^^   -^  372—   •  (4.4.39b) 

CO-x»  00 

The  plot  for  the  noise  spectrum  is  shown  in  the  1^=0  curve  of  fig.l, 

Spectrum  is  constant  at  low  frequencies  and  at  high  frequencies  the 

-3/2 
spectrum  goes  as  co    .  The  con^uter  program   used  to  plot  the 

spectrum  is  given  in  the  Appendix.  (u=wt  ,  ri=aT ,) . 


4.4.2  Diffusion  With  Sinks   (a^O) 

When  sinks  are  present  a  is  not  equal  to  zero.  The  required 

noise  spectrum  expression, taking  real  part  of  (4.4.37)  becomes 

16Tra^PQ(p^-5^)   letra^p^K 

S  p(M)   =  2   2  2   *   2~"TT  (4.4.40) 

^^         3D(p  +5^)      D(p  +6-  )^ 

where 


K  = 


and 


(10p^5^-5p'*6-5^)[2p5+25e'^cos5(p+2)+e"Psin5(5^-p^-4-4p)] 

- (p^-lOp  V+Sps"^)  [p^-6^i-25e"Psin5 (p+2)+e"Pcos6 (p^-5^+4-4p) -4] 

(4.4.41) 

„  ^  2a  [V^   +  ^Ja^    J  (4.4.42) 

Vd   ^____ 

g  _   2a  [V-^  +  jv'a^+co^  ].  (4.4.43) 
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Fig.  1 
Noise  spectra  for  embedded  sphere  (r|=0,  r|^0) 


51 


Furthermore  using  the  following  dimensionless  quantities 


2 
Tj  =  a  /D,     diffusion  time 


U  =       COT^ 


and 


(4.4.44) 


n        =      T^/T^      =         aT^ 


we  get  for  p  and  6 

p     =     V^"  ^/Tl+  Vt)  +u  (4. 4. 45a) 

and 

5     =     \^/-Ti+  VV]  +u         .  (4.4.45b) 

If  in  the  above  expression     r]  =  0,  we  have     p  =  5  =  y       and 

equation   (4.4.40)   reduces  to  equation   (4.4.38).     The  general  case 

a  ^  0     has  been  considered  for  the  first  time  as  far  as  we  know.     The 

low  and  high  frequency  asymptotes  are  found  to  be  the  same  as  for 

a  =  0 

32Trp^a 

^AP^"'^     "         15D  (4.4.46) 

co-K) 

and 

PQ(2D)^'^4Tra2 

^AP^"^     "    3/2 •  (4.4.47) 

a)-><»  CO 

The  noise  spectrum  has  been  plotted  in  fig.l .     For  low  fre- 
quencies the  spectrum  is   constant  and  at  high  frequencies,   the  spectrum 

-3/2 
goes  as  CO         .     For    ri  »  1     if  a  range  can  be  foiond  such  that 

3  2 

Ti  5  u  <  rj  ,  then  the  spectinm  goes  as  co       and  at  higher  frequencies 

-3/2 
CO  behavior  takes  over.     The  computer  program      used  to  plot  the 

spectrum  is  given  in  the  Appendix. 


CHAPTER  5 
N0NSY>*1ETRICAL  EMBEDDED  BODIES 


5 . 1   Introduction 


t»t.*->,. 


We  have  seen  in  chapter  2  that  for  nonsynunetrical  geometries 
Green's  function  cannot  be  found  in  closed  form  and  therefore  eigen- 
function  expansion  methods  should  be  used.      In  this  chapter  we  shall 
study  the  noise  spectra  arising  from  nonsymmetrical  embedded  bodies 
for  both  Voss-Clarke's        'P'   source  and  the  physical  diffusion  noise^^ 
source. 


5.2  Voss  and  Clarke's  'P'  Source 

5.2.1  Nonsymmetrical  Bar 

We  consider  the  diffusion  in  an  infinite  v-dimensional 
domain  with  fictitious  boundary  conditions  i.e.  embedded  case.  In 
all  integrations  of  chapter  2  d  r  is  replaced  by  d  r.  The  norm- 
alized eigenf unctions  of  (2.2.2)  are 

with  ^  =  (k^ , . . .  ,k  )  and 

2Trn. 
\     =  -^   ,  (5.2.2) 

i 

V 

n.  =  0,±1,±2,...  and  V  =  n  L.  ;  (5.2.3) 

^  i=l  ^ 

for  embedded  case  V  goes  eventually  to  «>.  The  eigenvalues  are 

2  2    2  2      2 

\^     =  M,  =  Dk  ,  where  k  =  k  +k  +...+k   .         (5.2.4) 
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Since  k-values  are  dense,  we  can  write 

2~    /d\  ZCy,     with     ZQQ   =  -y—  (5.2.5) 

h  (2TT) 

being  the  density  of  states  for  a  domain  of  dimension  v.  Since  all 

noise  expressions  are  bilinear  in  cp  or  \|*',  V  cancels  out  for  any  sum  Z 

k 
in  the  limit  V  -»■  ».  Then,  equation  (5.2.5)  is  exact  and  the 

eigenfunction  expansions  siirply  represent  Fourier  integrals.  Thus, 

(2.2.31)  results  in,  with  /  d  k  denoting  integration  over  infinite 

V  dimensional  ^- space. 


„bulk,  , 


^     1 1  d\  d\' 


C2TT) 


/     e^^--^  d^l    2 


(Dk^+iw)(Dk'^-ico) 


iV_       ,v 


X  //  d^r^  d-'r^  S^CC^,;^)   e 


■iQs--5i-is'X2) 


(5.2.6) 


For  Voss  and  Clarke's  P  source 


So  for  (5.2.6)  we  now  obtain 


(5.2.7) 


„bulk,  .          0         -,v, 
Sa        C")   = /d  k 

'^  (2Tr)^  « 


.       ik.r    ,v 
/     e  ~  ''^  d  r 

V 
s 

n2i,4     2 
D  k  +C0 


(5.2.8) 
L. 


Now  considering  V     to  be  a  rectangular  bar  of  dimensions--^  5  x.   <  ^r— 
i=l,2,..,,v,      the  result  for     S.  becomes 


^bulk,  ,  ^2,,,  ,2    r     d\         ^     """   ^-^ih^ 

h       f")     =     Pq^^s)     ^:2;T-2    .",      .1,    ,    .2 


(5.2.9) 


00  D  k  +00     i=l      (dc.L.) 
2  1  1-' 
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For  a  three  dimensional  rectangular  bar  v  =  3,  and 

2  1 
,  -.         -    -    ,3,    3  sin  C^.L.) 
„bulk.  ,     n2,,,  ,2-   d  k    „      2  i  i'  ,-   ^   ,^. 

h    f"^  =  ^o^^s^  ^  ;tt-t  .",    ,h  ,  ,2    •  f^-2-io) 

«  D  k  +m  1=1   f-^.L.) 
^2  1  i-^ 

For  a  three  dimensional  rectangular  bar  the  spectrum  shows  breakpoints 

2 
at  CO.  =  D/L.  ,  and  for  the  case  L. »  L^  »  L_  (nonsjrmmetrical) 

the  following  frequency  ranges  are  claimed  by  Voss  and  Clarke  , 

Qr   ^  -1/2 

CO   «    CO,  ,  b  (.COj      a       CO 


co^    «  00  «  co^,  S(co)    a  CO 

-3/2 

C0_    «    CO  «    C0_ ,  S  (co)      a:    CO 

CO       »    C0_,  SCco)      a:    CO  .  (5.2.11) 


Voss  and  Clarke   noted  in  particular  the  1/co  range  between  limits 
co^  and  co^ . 

In  an  effort  to  find  the  breakpoints  and  the  spectral  frequency 
ranges  the  method  of  numerical  integration  was  used.  The  complete 
analysis  is  given  below  and  the  listing  of  the  computer  program  which  uses 
the  trapezoidal  rule  for  numerical  integration  is  given  in  the  Appendix. 

Explicitly,  for  the  three  dimensional  bar  we  have 

,    ,,  00       dk     dk     dk 

S^'^'^'Cco)   =  const.   /  /  /       ^"^  ^  ^^     I  2     2  (5.2.12) 

-"       D    (k  +k  ,+k  )    +C0 
X     y     z-^ 

sin^^^ip        si"^(^yL2)        sin^C^^Lg) 
^Kh^  ^f^yS^  ^¥ih^ 
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Let  us  call 


V2     -     ^ 


k   L,     =     ^  (5.2.13) 

Using  the  above  transformation,  we  get  as   integrand 

dg  dr)  d^      1 

L1XL2XL3        ^  D^2  ^      D     2   ^     D^2^2^^2 

sin  -^rg  ^         sin  ^  -         sin  -^  _ 
X   C-y^)^  X   C-y^)     X   C-y^r    •  (5.2.14) 

2^  2^  2^ 


Furthermore  call 


coj   =   D/lJ,     CO2  =   D/L^,     coj  =   D/Lj  . 


u  =  co/Vw.'jO  00   , 


"1  =  ^"l/'^2"3  ' 


^2   =  ^"2/"l"3  > 


=  X^ 


(5.2.15) 


(5.2.16) 


"3  =  ^"3/"l''2  . 

We  then  obtain  the  integral  in  the  following  dimensionless  form, 
with  F(co)  =  const.  S(w): 


W- 


Fig.  2 

Noise  spectra  for  nonsymmetrical  bar 
(Voss  and  Clarke's  'P'  noise  source) 
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+  40 


♦  30  - 
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P,-,    ,  :,              dg  dri  iC                           '^"K2    ,^^"^2   ,^^"^2 
FCco)  =  /  /  /   ,   7  ■  '2  -  ^2,2  -2    <^-T ^^     ^1 ^  ^- ^ 


CWjC  *<^2^   +^3^  )  +C0      ^        ^        ^ 


(5.2.17) 


Finally  we  choose  as  parameters  L^  =  lO^m,  L„  =  1cm  and  L_  =  0.1cm, 

2  2 

D  =  25cm  /sec  (corresponding  to  mobility  of  1000  cm  /volt-sec.) 

Then  we  have  coj^  =  0.25rad,  co  =  25rad,  co  =  2500rad;   and 

-      -2   -  -      2 

coj  =  10  ,  00  =  1,  and  w  =  10  .  These  are  the  normalized  break- 
point frequencies.   (Units  for  oi's  should  read  as  Rad/sec.) 

Substituting  all  the  above  constants  and  using  the  symmetry 
property  of  the  integrand  we  have  for  (5.2.17) 

,_U      .1       .1 


p.--,  -   o    r   r   r         dg  dri  d^ r_l_^2   /"'T',2   ,1^,2 

0  0  0  (.oir+-n +ioor)  +^   4c      ^      k 


sin:^^  ^    sin^^  2    sin^ 
1  ^  ^~T 

(5.2.18) 


The  plot  of  the  spectrum  is  given  in  fig.  2.   We  have  not 

been  able  to  verify  results  as  in  (5.2.11).  The  spectrum  has  no  1/f 

range  and  goes  much  steeper  at  low  frequencies;  only  the  break- 
points at  CO  and  00  are  visible. 

5.3  Physical  Diffusion  Source 

We  have  seen  in  chapter  3  that  the  physical  diffusion  noise 
source  has  a  spectrum 

S^  =  4D  VV'nCc)6(x-X')-  (5.3.1) 

When  this  is  substituted  in  equation  (5,2.6),  which  is  the  general 
result  for  nonsymmetrical  embedded  diffusion,  one  finds  for  the 
spectrum 


iKf".-- 
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"■  2  1 

^A   Ceo)  =  C  /  -—-^    U  (5.3.2) 

~  D  k  +C0  1=1  Ct*-L.) 
2  1  1 

2 
which  differs  from  equation  C5.2.9)  by  the  factor  Dk  in  the 

numerator;  C  is  a  proportionality  constant.   (For  particle  diffusion 

C  =  4nQ/(2n)\) 


5.3.1     Rectangular  Bar 

For  a  three  dimensional  rectangular  bar  equation   (5.3.2)  becomes 

2   1 
u   11                       j3,    _,,  2     3     sin   (-dc.L.) 
„bulk,  .        r   r  ^  k  Dk       „                2  i  i'  ,-   _  _. 

S.        (w)   =  C  /  -y-^ — 2     ^     n^ 2 •  (5.3.3) 

"  D  k  +00     i=l    (dc.L. ) 
^2  1  1"^ 

Again  for  the  case  L^  »  L^  »  L_   (nonsjnnmetrical)  Voss 

and  Clarke  claimed  the  spectral  ranges  as 

f  «  f . ,  S(co)  =  constant 

fj  «  f  «  f^,  S(co)  cc  -log  f 

f2  «  f  «  f3,  S(co)  «  f'-^/^ 

f  »  fj,  S(co)  <x  f"^/^                      (5.3.4) 

Using  the  same  transformation  as  in  section  5.2,  the  additional 

2 
factor  Dk   is  found  to  be 

2      2      2 
=     D(k  +k  +k  ) 
X     y     z 

■•l     ^     H 

1/3  -      2   -      2   -      2 
=      (cOjCo^Wj)        (co^?  +C02T1  +00^^  )    .  (5.3.5) 
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The  overall  noise  spectrum  is  given  as 


oo  oo  oo 


F(^)   =   C  /  /  /  5 '-^ %   „  % 

0  0  0  e^^r+m^^  +^34  )  +^'' 

.1^        .1        .    1^ 
smy?      sm^      sm^ 

X  (-^^-)   X  i-^^y   ^  (-j-^)    .  (5.3.6) 

2^  2^  ^ 

2 
For  the  values  of  L^  =  10cm,  L„  =  1cm,  L_  =  0.1cm,  and  D  =  25cm  /sec. 


this  becomes 


F(«)  =  Cf  f  f 


dKdr]dK   (.OlC^+ri^+100^^) 


2  2     2  2-2 
(.Oir+ri^+100^  r+M 

.  1„       .1        .    1^ 
sm^      smyH  j         sm^ 

x(^-^)  X  (-pA-)^  X  i-^)^    .  (5.3.6a) 

2?        2^        2^ 


The  listing  of  the  computer  program  is  given  in  the  Appendix. 
We  have  not  been  able  to  verify  the  results  as  claimed  by  Voss  and 
Clarke's  equation  (5.3.4).  The  spectrum  is  given  in  fig.  3. 
The  high  frequency  slope  is  -3/2  as  predicted  by  Lax  and  Mengert  • 
Breakpoints  oj  and  oo_  are  visible  and  at  low  frequency  (not  shown) 
the  spectrum  levels  off.  Also  there  is  no  obvious  1/f  range  in 
the  spectrum. 

Furthermore  we  can  also  see  that  in  the  noise  spectrum  for 
the  nonsymmetrical  bar  for  both  the  noise  sources  (namely  'P'  source 
and  diffusion  source),  figs. 2  and  3,  at  the  low  frequencies  there  is 
a  "wiggle"  in  the  spectrum.  The  presence  of  such  a  "wiggle"  must  be 
attributed  to  the  nonsymmetry  of  the  volume  under  consideration. 


-o;i 
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Noise  spectra  for  nonsymmetrical  bar 
CPhysical  diffusion  noise     source) 
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The  "wiggle"  appears  on  the  noise  spectrum  towards  the  longest  dimension 
of  the  specimen.  As  we  shall  see  in  the  next  subsection,  for  more  sym- 
metrical geometries  like  square  bar  and  cube,  the  noise  spectrum  is  smooth 
and  the  low  frequency  "wiggle"  disappears  (figs.  4  and  5  ). 

5.3.2  Square  Disc 

The  noise  spectrum  arising  from  a  square  disc  has  been  analyzed. 

In  equation  C5.3.6)  the  parameters  have  been  changed.  The  longest  dimen- 

2 
sion  has  been  eliminated  and  a  square  disc  of  1cm  cross  section  and 

thickness  of  0.032cm  has  been  considered.  The  noise  spectrum  is  given  as 
0    (O.ir+O.lTi  +100^  )+w 


.1^       .1        .  L, 
sin-H  ^    sin-rT)     sm^  - 

^{.-Y^y  ^  (-1"^-)  X  i-Y^y        •  (5.3.7a) 

^  2^  -f 


The  computer  program  listing  is  given  in  the  Appendix.  The 
plot  for  the  spectrum  is  shown  in  fig.  4.  We  can  see  that  the  spectrum 

varies  smoothly  and  the  low  frequency  "wiggle"  disappears.   Further 

-  3/2  -   ~ 

at  high  frequency  the  spectrum  goes  as  f    .  Breakpoints  atWj^=W2  and 

u,  are  visible  and  at  low  frequencies  (not  shown)  the  spectrum 

levels  off. 


■■««! 


Fig.  4 

Noise  spectra  for  square  disc  and  square  bar 
CPhysical  diffusion  noise  source) 
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5.3.3     Square  Rectangular  Bar 

The  noise  spectrum  resulting  from  a  square  rectangular  embedded 

bar  has  been  studied.     The  smallest  dimension  has  been  eliminated 

2 
and  a  square  bar  of  crosssection  0.1cm     and  length  3.16cm  has  been 

considered.     Equation   (5.3.6)    for  the  noise  spectrum  is 
0        (O.Oir+lOri^+lO^  )+co 


.    1„  .1  .1^ 

sin-=^  sm^r^  sin^ 

x(-T^)  (.-Y-r  (.-^r    .                                                      (5.3.7) 

2^  fl  2^ 


The  computer  program  listing  is  given  in  the  Appendix.  The  plot 
for  the  spectrum  is  shown  in  fig. 4  .  We  see  that  the  spectrum  varies 

smoothly  and  again  the  low  frequency  "wiggle"  disappears.  At  high 

-3/2  -  - 

frequencies  the  spectrum  goes  as  m   .  Breakpoints  at  w,  =03.  and  co 

are  visible  and  at  low  frequencies  (not  shown)  the  spectrum  levels 
off. 


5.3.4  Cube 

The  noise  spectra  obtainable  from  a  less  asymmetrical  embedded 
geometry  like  a  cube  has  been  studied.  The  longest  and  the  smallest 

dimensions  of  the  nonsymmetrical  bar  have  been  eliminated  and  we 

3 
considered  a  cube  of  dimension  1cm  .   The  noise  spectra  is  given  by 

sr^^)  -osff  ff¥4^  c|ii,^  c^,^  ^^  . 

-      (.K  +11  +^  )    +w  -C  j'H  2^ 

(5.3.8) 
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The  plot  for  the  noise  spectrum  is  shown  in  fig.  5.  We  can 
see  that  the  spectrum  is  very  smooth,  the  spectrum  levels  off  at 
low  frequencies  and  the  high  frequency  slope  is  f  '  .  The  break- 
point at  co_  is  quite  visible  and  the  roll  off  is  much  steeper  than 
for  the  square  bar.  The  computer  program   listing  is  given  in  the 
Appendix. 

Thus  in  this  chapter  we  have  analyzed  the  noise  spectra 
obtainable  from  nonsymmetrical  embedded  bodies  for  'P'  source  and 
nonsymmetrical  as  well  as  less  asymmetrical  embedded  bodies  for 
the  diffusion  noise  source.  We  note  that  the  spectral  ranges  claimed 
by  Voss  and  Clarke  are  not  verified.  For  the  'P'  source  their  claim 
for  a  1/f  noise  spectrum  is  far  from  true.  Further  we  have  seen  that 

the  noise  spectra  for  the  physical  diffusion  noise  source  results 

-3/2 
in  f     type  of  behavior  at  high  frequencies.  There  may  be  a  quasi 

frequency  range  where  1/f  kind  of  behavior  is  observed,  but  this 
range  is  small. 

In  the  next  chapters  we  shall  study  the  effect  of  the  boundary 
conditions  on  the  noise  spectra,  both  for  the  symmetrical  and  non- 
symmetrical unembedded  geometries. 
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Fig.  5 

(Slope  indicated  is  per  decade  of  frequency) 

Noise  spectra  for  cube 
(Physical  diffusion  noise  source) 


CHAPTER  6 
SYMMETRICAL  NONEMBEDDED  BODIES 


6.1   Introduction 


In  this  chapter  we  shall  study  the  noise  spectra  resulting 
from  symmetrical  nonembedded  bodies.  We  already  have  discussed  in 
chapter  1  the  difference  between  embedded  and  nonembedded  cases, 
the  latter  being  subject  to  boundary  conditions  at  the  surface  of 
the  body.  We  thus  investigate  the  effects  of  the  boundary  condi- 
tions on  the  noise  spectra  for  symmetrical  geometries.  For  example, 
in  the  case  of  surface  generation-recombination,  we  have  a  mixed 
Dirichlet-Neiimann  boundary  condition  for  the  fluctuating  carrier 
density  An, 

aAn  +  bVCAn)  =  0  C6.1.1) 

where  a  =  s,  surface  recombination  velocity  and  b  =  D,  diffusion 
constant.  The  presence  of  boundary  conditions  results  in  (i)  altering 
the  shape  of  the  spectra  from  that  of  the  "infinite  medium"  and 
(ii)  correlation  effects  may  occur. 

6.2  Summary  of  Previous  Work 

In  this  section  we  briefly  review  the  results  of  the  work  done 

so  far  to  obtain  noise  spectra  from  symmetrical  nonembedded  bodies 

such  as  linear  bar  and  cylinder. 
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6.2.1  Linear  Bar 

The  problem  of  one  dimensional  diffusion  in  a  finite  domain 
was  studied  by  van  Vliet  and  Fassett^.  The  one  dimensional  diffusion 
in  the  region  -L  5  x  5  L  with  perfectly  reflecting  barriers  at 
-L  and  +L  is  considered,  where  L  >  d.  The  spectral  density  for 
the  fluctuating  variable  AP,  S  C")  in  the  subregion  -d  5  x  <  d 
is  foimd  to  be 

3 

S      (0,)    -  ^0^^^^      sine-Hsinhe-FCa9.be) 

^AP*-  -•      D      3  (6.2.1) 

e-'ccoshe-cose) 

where  p   is  the  average  number  of  particles  per  unit  length, 

F(ae,be)  =  sinae  coshbe  +  cosae  sinhbe  +  sinbe  coshae  +  cosbe  sinhae 
and 

e  =  2L\/w/2D,    a  =  p   b  =  ^  . 

Equation  (6.2.1)  reveals  that  at  sufficiently  high  frequencies  the 

-3/2 
spectrum  varies  as  w   ,  whereas  at  low  frequencies  the  spectrum 

becomes  constant 

—  (2D) 

^AP^'^^   =   ,,,5  3/2  (6.2.2a) 

o)-**   (2L)  CO 

and 

S^pC"^   =       '/        ■  (6.2.2b) 

a>^     3L 

Thus  the  presence  of  reflecting  boundaries  limits  the  low 
frequency  noise,  whereas  without  walls,  as  in  the  embedded  case  the  spec- 
trum rises  indefinitely  as  co  ->■  0  according  to  co  '^. 
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6.2.2  Cylinder 

The  case  of  f luctuatiore  in  the  finite  cylindrical  domain  with 
mixed  Dirich let -Neumann  conditions  was  investigated  by  van  Vliet- 
Chenette^^  in  1965.   Following  Morse  and  Feshbach^**,  it  was  shown 
that  the  Green's  function  or  consequently  the  total  spectra  can  be 
written  as 

gtotal  ^  s"  -.  S^  (6.2.3) 

where  S  ,  is  the  spectrum  as  in  the  embedded  case  and  S  is  due  to 

the  presence  of  boundary  conditions.  Following  this  procedure  the 

total  spectra  for  the  nonembedded  cylinder  was  found  to  be 

2 
gtotal  ^  s°°  +  S^  =  4(varX)^  -Re  P(co)  (6.2.4) 


where 


P(j^)     = ^{1  +  — 


,  ,2     xa  DxJ,  (xa)-aJ-(xa) 
(xa)  1      0^ 


J's  are  the  Bessel  functions 

a  is  the  radius 


ivi^ 


D 

and  a  is  a  constant; in  the  case  of  surface  recombination,  a  is  the 
surface  recombination  velocity. 

The  spectrum  is  constant  at  low  frequencies  and  at  high  fre- 

-2 
quencies  the  spectrum  goes  as  oo~  ,  with  the  presence  of  a  diffusion 

-3/2 
range  (oo    )  depending  on  related  time  constants.  Thus  due  to  the 

presence  of  boundary  conditions  the  logarithmic  low  frequency  behavior, 

as  in  the  embedded  case  disappears  and  the  spectrum  becomes  flat. 


Ii*5^^.:  .  .■'  :,.'.^-.'^^ 
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6.3   Present  Work:  Sphere 
In  this  section  we  shall  consider  the  case  that  the  domain 
V     =  V     is  a  spherical  domain  of  radius  'a'  with  homogeneous  boundary 
conditions  on  the  boundary  of  the  domain 

d|7G(x^,s,^  )  +  aG(x^s,X  )  =  0  (6.3.1) 

where  r,  is  a  surface  position  coordinate,  n  is  the  outer  normal 
and  a  is  the  surface  recombination  velocity.  The  differential  equa- 
tion for  the  Green's  function  is 

d^G       2   dG          2^^    1,   1        3  ,    .    o  d   .  1        9^   t^  1^, 

— T  +  -  J—  +  ■^t  G+-T-r— r— 5"  -r^  (sine  -tk)   +  ^ ^]G  =     -  Tr^fr-r') , 

J   2        r  dr                  2'-sin9   39  dB-^  .    2„   ^   2-'  q    v^ -</  -- > 

dr                                r  sin  9  Sep 

(6.3.2) 
where     G   C^,^' )    =  G(^,ioo,r'),   and  x  is  the  same  as  in  equation   (6.2.4) 
The  polar  axis  we  take  along  ;r' . 

Following  the  treatment  as  in  Methods  of  Theoretical  Physios, 
vol.1,  by  Morse  and  Feshbach      ,  page   825,  we  expand  the  Green's 
function  in  terms  of  a  complete  set  of  functions  involving  all  but 
one  of  the  coordinates    (9     and     ep)   with  coefficients  which  are  undeter- 
mined functions  of  the  uninvolved  coordinate.     Thus  let 

G^(X>X')      =       2     p^''"(r,r')p'J(cose)ei'"('P-'P')  (6.3.3) 

t     =  0,1,. . .,  »     ;       m  =  -i,..  .  ,+e  . 

Substituting  equation  (6.3.3)  in  (6.3.2)  for  the  first  three  terms  on 

the  left  hand  side,  we  get 

,2  t  ,m  J  i ,m 

V  r  '^       ^  2  ^x   ^2  i,m-,   m,   „-.  im(cp-(p')         ,.  _  .. 

Z[ 2~      T~d^     +  X  p^'  ]p^(cose)e   ^^  ^  ^         (6.3.4) 

^,m   dr 
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and  the  remaining  terms  on  the  left  hand  side  of  (6.3,2)  give 


25 


{Quantum  Medhanics ,  Schiff     ,  page  75) 


1    r     1     9    r-.inO  -^1   +        1       9     i   m,         .    imQcp-cp') 
2Uine  90 f-^^"®  W   *     .2^  — jJp^(cose)e 
r  sin  6  d(p 


— ^  '^  p    (cos9)e 


im(cp-(p*) 


2       ^m* 


(6.3.5) 


Now  expanding  the  delta  function  on  the  right  hand  side  of  equation 
(6.3. 2) ,  we  get 


srr- rn  -  5(r.-r')  5(cp-(p')5(e-9') 


sin9 


(6.3.6) 


Let 


imcp 


(6.3.7) 


-2  25 

where  (N.  )    is  the  normalization  constant;  according  to 

{Quccntum  Medhanias ,   Schiff,  page   73,    equation    (14,16)) 


N    =  /mnriE 

tm  4Tr   (^+ 


m 


m 


)! 


(6.3.8) 


Then  by  inversion  from  (6.3.7)  we  get 


F(e-,<p')  =  /Sie^:i|ip:iN^^p;(cose)e— dJ2 


•im(p 


(6.3.9) 


where  dS2  =  sine  d6d(p  ;thus 


FCe-.cp')  =  N^^pJ'(cose')e"^'"'<''  . 


Hence  we  get 


5U-^')=n2   2  p'n(eose)p';(cose')ei'"('P-<P')  ^-^^     (6.3.10) 


4 
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Since  the  polar  axis  is  along  r',  6'  =  0  and 

P'j(cose')  =  1.  (6.3.11) 

The  result  of  the  full  substitution  and  equating  term  by  term, 
yields 

,2   I ,m  ,   t ,m 

^x       ^  2  ^Px  2,       ^jl*!)  J,m  _  5(r-r') 

"T2~  *  tST-  *  ""   f^ r-JPx  -•      2„2        •  (6.3.12) 

dr  r  Dr  N. 

The  solutions  of  the  homogeneous  equation  corresponding  to  (6.3.12) 
are  given  by  spherical  Bessel  functions  of  first  and  second  kind 
j-(xr)  and  y»(xr)  respectively. 

The  complete  solution  of  equation  (6.3.12)  is  given  as 

''w^j(r)w^jj(r'),    r  5  r' 


""      r^D[W(w^j,w^jj)] 


V. 


(6.3.13) 
w^j(r')w^jj(r),   r  >  r' 


where  W  is  the  Wronskian  evaluated  at  r=r'  and  w   and  w    are 
two  independent  solutions  which  satisfy  the  boundary  conditions  at 

r=0  and  r=a,  respectively.  For  r=0  the  solution  must  be  regular.  We 

thus  have: 

w^j  =  j^(xr)  (6.3.14) 

and 

w^jj  =  Aj^(xr)  +  By^(xr).  (6.3.15) 

w.^^  satisfies  the  boundary  condition  (6.3.1)  which  is 

D(^)^  =  -o^.U^a).  (6.3.16) 
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Substituting  w    in  equation  (6.3.16)  we  get 


DA[ 


9j^(xr) 


8r   -'r=a 


]    +  D 


y^Cxr)^ 


ar 


r=a 


-a[Aj^(xa)+By^(xa)]    (6.3.17) 


or 


DB 


ay^ (xa) 


a j  ^ (xa) 
9 J.    +  CTBy^(xa)  =  -DA  — aAj^(xa) 


(6.3.18) 


Substituting  the  recurrence  relations,  see  ref.  26, 


ay^ (xa) 


ar 


y^  C'ta) 


and 


aj^(xa)     ^ 


ar 


a  ^^'^^^^ 


xy^^j(xa) 


xj^,^(xa) 


(6.3.19) 


(6.3.20) 


in  equation  (6.3.18)  and  simplifying,  we  obtain 


B 
A 


D[-  j^(xa)-xj^^^(xa)]  +aj^(xa) 


D[-  y^(xa)-xy^^^(xa)]  +ay^(xa) 
A  suitable  choice  for  the  coefficients  A  and  B  is: 


(6.3.21) 


m 


A  =  if(a+  — )y^(xa)    -  Dxy^^j(xa)J 


and 


D^,  . 


B  =   -i{(a+  f^)j^(xa)    -   Dxj^^j(xa)]. 


(6.3.22) 


(6.3.23) 


Substituting  A  and  B  in  equation    (6.3.15)   we  obtain  for  w..^   as 

^11     =     j^(xr)r(a+— )iy^(xa)-Dxiy^^j(xa)] 
-  iy ^  (xr)  [  (cj+— )  j  ^  (xa)  -  Dxj  ^^  ^  (xa)  ]    . 


(6.3.24) 
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Now, 


-y^Cxa)  =  h^  ■'(xa)  -  j^(xa) 


(6.3.25) 


(1) 


where  h^  (xa)  is  the  .«?pherical  Hankel  function  of  the  first  kind. 
Thus  we  obtain 


^^11  =  J^(^rH(a+  ^)h^^^(xa)-Dxh|Pj(xa)]-h|^^(xr)[(a+  ^)  j^(xa)-Dxj^^j  (xa)] 


which  can  also  be  written  as 


w^jj  =  h^  -"(xr)  -  j^(xr)Q^(xa) 


where 


Q^Cxa)  = 


n  u(l)^  ,  ,  D^,,  (1)  ,  . 


Dxj  ^^  ^  (xa)  -  (cj+— )  j  ^  (xa) 
The  Wronskian  as  defined  by 


(6.3.26) 


(6.3.27) 


(6.3.28) 


W  = 


'^I    ^^I ' 


^^11   w^II' 


(•  denotes  first  derivative)  evaluated  at  r  =  r'  is  found  to  be 


"f"^r««ii^ 


xr 


,2 


(6.3.29) 


Substituting  for  r  <  r'   in  equation  (6.3.13)  we  get 
p;'"(r,r')  =  — ^  [j^(xr)h^JJ(xr')  -  j^(xr)  j^  (xr')Q^(xa)l  . 


(6.3.30) 


Finally,   substituting  equation   (6.3.30)   and  equation   (6.3.8)   in 
equation    (6.3.2)   we  find  for  the  Green's  function,    for     r  5  r' 
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X  p'J(cose)e^'"^'^-'P'^  .  (6.3.31) 

For  the  noise  spectrum  we  must  evaluate  the  six  fold  integral 

S(co)  =  2Re/---/   r^dr  dCcose)dcpr'^dr'd(cose')dip'G  (x,;C')    (6.3.32) 
6  fold 
reduced 
domain 

where  the  integral  involves  the  reduced  domain  0  <  r  <  r* ,  0  S  r'  <  a. 

Integration  over  cp  gives  2Tr  for  m=0;  it  gives  zero  for  m?^0.  Next  we 

use  the  addition  theorem,  ref.25  which  gives  the  relation 

^j^    =   Z     (2^+l)j^(xr)h^^^(xr)P^(cos0).  (6.3.33) 

Substituting  in  equation  (6.3.31)  and  separating  the  noise  spectrum 
into  the  two  parts,  indicated  by  the  [  ]  of  6.3.31,  we  find  for  the 
first  part 

ixR 

S    (w)=       4tt  Re  /•••/     r  dr(dcosG)d(p'    r'    dr'd(cose')  j-^  (6.3.34) 

^  5   fold 


and 


„    ,  .         4var  X  =    ,  . 
Sj(w)     =  — y Sj(co). 

s 


The  equation  (6.3.34)  for  S  (oo)  is  exactly  the  same  as  for  the  spectrum 
of  the  embedded  sphere  computed  in  section  4.3;  this  spectrum  will 


now  be  denoted  as  S„(cjo) 
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The  second  part  of  the  spectrum  has  an  integrand  which  is 
regular  for  r  =  r' ,  the  solution  being  the  same  for  r  5  r'   as 
r  >  r'.  Hence,  we  can  at  once  integrate  over  the  full  domain 
0  S  r  S  a,  0  5  r'  5  a.  Thus, 


SjjCc^)     =  -Re  ^^/  •••  /  r^drdCcose)dCcose')r'^dr'd(p' 

5  fold, 
full 
domain 


X     Z      (2^+l)j^(xr)j^Cxr')Q^Cxa)P^(co5e).  (6.3.35) 

^=0 


The  integration  over     d(cos6')dcp'      gives   -4n.      Further  due  to  the 
orthogonal  property  of  Legendre   functions,   integration  over  d(cos6) 
contributes   -2  only  for  ^=0: 


/  P^(cose)d(cose)   =   -/  P^(z)PQ(z)dz  =  -25^p. 

Thus  we  have 

SjjCco)   =-Re-ip^//  r^dr  r'^dr'   j^Cxr) JQ(xr')QQ(xa)  (6.3.36) 

(Q^(xa),   substituting  t=0  in  equation    (6.3.28).     We  can  separate 
the  integrations  over  r  and  r' ;   thus  evaluating  the  integral 

frdTJ^dxT)      =     fr^dr  -^^ \/|     .  (6.3.37) 

Vxr 

In  equation  (6.3.37)  we  have  used  the  relation  which  connects  regular 

2  6 

and  spherical   Bessel   functions  J  and  j,   which  is      > 

\f^^     =  ^    Jn.1/2    f^^    •  ^6-3.38) 


78 


We  can  write  equation    (6. 3.. '^7)    as 

/r2drj^(xr>  r  2L£!lC-l2£lZ2ix^    ^    ^  (6.3.39) 

0  v'xr  X 


^^4    {      (xr)^/2dCxr)J^/2'^xr)  (6.3.40) 


xa 

/ 
x"     "     0 

We  can  also  use  the  derivative  formula  for  the  Bessel  function  which  is 


^[z\(z)]  =  z^J^.iCz)  (6.3.41) 

this  gives  for  v  =  3/2  the  required  integral  as 

/  z^f^J^^^{z)dz   =  z^/^J2/2(z)  (6.3.42) 

which  finally  gives  for  the  integral 

^^^^^^^0^^^^=  7-7372  4^2^^^^    .  (6.3.43) 

(xa) 


And  similarly  we  get  for 


3 


/r'2dr'Jo(xr-)  =  ^^  ^|  j   (xa) .  (6.3.44) 

(xa) 

Substituting  equations (6. 3. 43)  and  (6.3.44)  in  equation  (6.3.36) 

we  get 

.  .  3 
Sjj(co)=-Re  -^^  |[J3/2(xa)]'^QQ(xa)  (6.3.45) 

and  for  the  spectrum  we  obtain 

r.     r   -N     4varX  =   ,  . 
s 


->■<«( 
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4     3 
Taking     V     =  -^a  ,     we  obtain  for  the  fluctuating  variable  X, 

c      /■  1  n       6TrivarX   ft       ,•     ^^2^    ,     . 

SjjCco)  =  -Re  2—  [Jy2^^^^^  %ix-S^)  ■  (6.3.46) 

ux 


can  write  the  total  spectrum  for  nonembedded  sphere  as 


Denoting  the  spectrum  associated  with  second  term  by  S^(co),  we 


-total-  ,     _«.  .   _A,  - 

X*^"^  "  ^X^"^  ■"  ^x^"^'  C6.3.47) 

Now  we  shall  derive  another  form  for  s"(u)  in  terms  of  Bessel 

A 

functions  more  suitable  for  combining  with  S  (co) .  From  equation  (6.3.34) 
we  can  write 

Sj  =  Re  — /         /     r  dr  r'   dr'   d(cose)d(cose')   dcpdcp' 

s  6  fold, 
reduced 
domain 

X     2      -i^j^Cxr)hJl^(xr')p'J(cose)e^'"'^^-'P'^      ;  (6.3.48) 

^m     DN^^ 

integration  over  dcp  gives  2n   for  m=0  only, 
integration  over  d(cos6)  gives  -2  for  1=0   only  and 
integration  over  d(cos9')dcp*  gives  -4Tr. 
Thus 

=  Re  l^'f  ^;  r'^dr.h^'cxr')  J'Arj.Cxr).      C6.3.49) 
3  na    ^   0  0 

These  integrals  are  similar  as  found  in  equation  (6.3.37),  thus 


r'   2  ,3/2     - 

/rdrJQ(xr)     =     I_- Vf  J3/2  (xr-) .  (6.3.50) 


W: 
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So  we  get  for  S^ 

S       -     r.-   8varX  4TTix  rr  ^   (xr')    dCxr')Cxr']     nCD  c^^n j        fvrn      f6  3  511 
b        -      Re  -^^  D^2i        6,      ,,3/2/— r         "1/2'-^''   ^•^3/2^'^''  ^       C6.3.51J 

s  Q     X   Cxr']   '   vxr'  '  ' 

p 

=  Re  — r^ D~~5l-^  ^^^^  d(xr)HJ-,^(xr)J  ,  (xr).         (6.3.52) 

s     X    0  ' 

2  7 

The  required  integral  was  calculated  by  Luke 

a  4 

/  (xr  )^dCxr  )Hp^Cxr  DJ^Cxr  )  =  -^^^{h[^^  (xa)  J^Cxa)  +  J^Cxa^H^^^  (xa)}- 

°  11  2"      2"       F    2 

And  using^^  (equation  9.1.27)  which  is  (6.3.53) 

J5/2Cxa)  =  -Ji/2Cxa)  +^3/2Cxa),  C6.3.54) 

and  further  using  the  Wronskian  relation  we  foiind  that 

'^3/2^^^)"}/^'^^)  -  ^1/2^4/2^^^)  =  HI  •  f^-3-55) 

Substituting  equations  (6.3.54)  and  (6.3.55)  we  found  that  for 
equation  (6.3.53)  (also  using  (6.3.38)): 

=  -^^Sl^f  ^l(-^<^(-^)   •  (6.3.56) 

For  S^  we  obtain  finally 

Sj(co)  =  Re^  f-  fci"*  i  jj(xa)hP^(xa)]  (6.3.57) 

where 


fe»^'-.^- 
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The  spectrum  for  S^^    from  equation   (6.3.46)    is 

Sjj  =-Re  ^ [J^^^ixa)]   QqCxe)  (6.3.59) 

ux 

„   12fvarX)a   r,/^     ,        ,     ,,2-    ,     , 


"^^f^l'^''^^^\'^''^^  (6.3.60) 


where  QQ(xa)  from  equation  (6.3.27)  for  ^  =  0 

DxhP^(xa)-ah^^^(xa) 
^O^'^^^  =  Dxj^(xa)-aJQ(xa) '  ^^'^-^D 

For  the  total  spectrum  we  now  find 
total       A     1  rn       7     ^^1   (>ta)-ah^  '' (xa) 


(6.3.62) 


Solving  for  {    }  quantity  first  we  get 


-CThJ  Uxa)JQ(xa)+ah|^^^(xa)jj(xa) 


'^    >     =     '^1^^^^     Dxj,(xa)-aJo(xa) '  C^'^'^^) 

Using  ref.  25  (equation  (10.1.31))  we  have 

ho^^(xa)j^(xa)  -  h^^  (xa)  JQ(xa)  =  — ^-^  .  (6.3.64) 

(xa) 

So  we  get  for  {  }  quantity 

J        CT j .  (xa) 

Thus  we  get  for  the  total  spectrum: 
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-total,  ,  ,n         1        rl       ^  JjC^a) 


Further  we  know  that 


xa 


=     iaV2^  (6.3.67) 


(xa)   =  -a  (-5—).  (6.3.68) 

Introducing  the  following  time  constants 

T   =  —  (volume  lifetime)  (6.3.69) 

2 
T.  =  —  (diffusion  time)  (6.3.70) 

T   =  —   (surface  lifetime)  ,  (6.3.71) 

and  also  introducing  the  dimensionless  quantities - 

U  =  COT^,    Tl  =  T^T^,     C  =  T^/Tj  (6.3.72) 

we  can  write 

2 
(xa)   =  -(iu+Ti).  (6.3.73) 

We  can  rewrite  equation  (6.3.66)  for  S^°  ^  as 

-total,  ,     Td   1   ri   3ct  h^^^^  ,  ,,   ,  ,,, 

^X    Ceo)   =  ARe^^[U-^^^^^^^_^^^p^^.^^^^^]^     (6.3.74) 

2 
where  A  =  —  =  — ^^ — - — ^ —  .    Using  the  time  constants  and  dimension- 
less  quantities  defined  in  equations  (6.3.69)  to  (6.3.72),  we  have 

5^  =  -  3==  '  (6.3.75) 

and 


1 

fyrj+iu 

1 

Dx  a 


2^   =  -  ?(Ti+iu)  •  (6.3.76) 
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Using  equations  (6.3.75)  and  C6.3.76),  the  spectrum  finally  becomes 

gtotal ^^^  ^  4(var  X)a^  ^^p  ^^^  ^^  3  _  ^^^ 


where 


PCco) 


^'   Ur] 


3  j  ^  CiVr)+iu)  1 

+  iu) j  (iVr)+iu3  +  iVin+iu  j'  CiVri+iu)  J  *  tD.3.78) 


■q+iu 
Now  we  shall  discuss  various   cases. 

6.3.1  No  Volume  Sinks 

If  there  are  no  volume  sinks  present,  then  r]  =  0.     For  this  case 

(■        3j  CiViH)  1 

RePCo))  =  Re  ^1  -  ^,.    ,  .  t-.f^^    •,/■ ,^WS   •      (6.3.79) 

lu     5(iu)j  (ivTu)  +  iviu  J  (iViuJ  I        ^      -" 

The  first  term  does  not  contribute.  Using  the  relations  for  spherical 
Bessel  functions    we  have 

jgCz)   =  ^  (6.3.80) 

and 

.    .  .            sinz       cosz  ,,    _   „,^ 

JjCz)     =     —2 —  •  (6.3.81) 

z 

Thus  we  obtain 


4a"?(C  +D  ) 


where   a  =  vy 


A  =  (-sina  cosha+a  cosa  cosha+a  sina  sinha)  (6.3.83) 

B  =  (cosa  sinha  +  asina  sinha^a  ;osa  cosha)  (6.3.84) 


/•;> 


^ 
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C  =  (cosa  sinha+a  sina  sinha-a  cosa  cosha  - )        (6.3.85) 

D  =  (-sina  cosha+a  cosa  cosha+a  sina  sinha  +  liIi°L_£2L_EL)  ,       (6.3.86) 


At  low  frequencies,  u  «  1,  and  the  spectrum  reaches  a  limiting 
value.  For  finite  domain  a  plateau  is  always  reached.  Expanding  the 
trigonometric  terms,  we  find  canceling  terms  up  to  order  a  : 

Re  P(co)  =  ll|i  .  (6.3.87) 

CCr+0 

At  high  frequencies,  u  »  1,  and  the  spectrum  reaches  the 
asymptote 


Re  P(co)  =  — -  .  (6.3.88) 

For  a  closer  investigation  we  must  distinguish  between  two  cases. 

6.3.1.  (a)    .     First,  let  C  «  1,  or  t  «  t..  This  option 

s     d        ^ 

corresponds  to  the  "diffusion  limited"  case.   In  equation  (6.3.87) 
the  terms  in  ?  are  now  irrelevant  and  the  spectrum  reaches  the  low 

frequency  limit  Re  P(0)  =  yr  •  At  high  frequencies  there  is  another 

-2  2 

region,  prior  to  the  u   region,  in  the  range  1  «  u  <  1/?  . 

6.3.1.  Cb)    .     The  surface  limited  case  occurs  when  ^  »  1, 
or  X     »  Tj.   In  this  case  the  diffusion  process  is  practically 
instantaneous  and  the  noise  is  only  due  to  surface  sinks.  At  high 
frequencies  the  limiting  behavior  is  —^  for  any  u  »  1.  However, 

this  spectrum  is  also  <^ound  at  small  u,  since  in  the  range 

1  3 

Y<   u«   1,     we  also  found  the  limiting  value  of  the  spectrvim  to  be  — =-. 

1  ^" 

This  biahavior  occurs  therefore  for  all  u  >  — .  The  turnover  occurs 
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for  u  «si  —  or  CO  sa  —  .  The  plateau  that  then  is  reached  has  the 
value  ■=■. 

The  full  spectrum  for  various  C  has  been  plotted  and  is  shown 
in  fig. 6.  We  see  for  ?  »  1,  or  for  the  surface  limited  case, the 
turnover  frequency  is  u  =  2/C.  For  ^  «  1,  or  for  the  diffusion 

limited  case,  for  K  =   .01,  the  diffusion  type  behavior  is  observed, 

3 
in  the  range  10  <  u  <  7x10  .  The  computer  program  used  to  plot 

the  noise  spectrum  is  given  in  the  Appendix. 

6.3.2  No  Surface  Sinks 

If  CT  =  0,  —  =  0.  This  corresponds  to  S,   being  infinite. 
The  spectrum  despite  the  possible  presence  of  microscopic  diffusion 
exhibits  only  a  regular  relaxation  spectrum. 

Re  P(co)   =  -y3_  (6.3.89) 

Tl  +U 

or  the  spectrum 

„f  ^           4Cvar  X)t  ^,   ,  __., 

SCw)     = V     .  (6.3.90) 

1+00   T 
V 

This  is  the  "volume  limited  case". 


6.3.3  Both  Surface  and  Volume  Sinks  Present 

When  both  surface  and  volume  sinks  are  present  we  have 
that  "nj^O.  We  again  use  the  relations  for  j^  and  jr,  as  given  in  equa- 
tions (6.3.80)  and  (6.3.81)  and  after  much  algebra  the  real  part  is 
found  to  be 


'•^v..^ 
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10*        10'        10*        10* 


Fig.  6 


Noise  Spectra  for  Nonembedded  Sphere  Cri=0) 

-3/2  -2 

(Slopes  indicated  are  u      and     u  ) 
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where 


P     =     V^  +  j\/ri+u^  ,  C6.3.92) 


,rT 1  /T     2 

=       V-    ^Tl    +    -VVl    +U 


5     =     V-  yTi  +   ^-VTl  +u  ,  (6.3.93) 


A     =      C-sin5coshp+6cos5coshp+psin5sinhp)    ,  (6.3.94) 

B     =      (cos5sinhp-pcos5coshp+5sin5sinhp)    ,  (6.3.95) 

C     =     A,HIl5cosM       ^  (6.3.96) 

D     =     B  -  £25|5inhi       ,  (6.3.97;, 

P     =      (p^-5^)C-2p5D  (6.3.98) 

Q     =      (p^-5^)D+2p5C  (6.3.99) 

M     =      (P^-6^)P   -   2p5Q  (6.3.100) 

N     =      (p^-5^)Q+2p5P.  C6.3.101) 

At   low  frequencies  in  the  above  equation    (6.3.91)  p  =  i/n 

and     5  =  0.     The   low  frequency  plateau  value  for  small  r)  is   found 
to  be 

Re  P(co)      =     i^  (6.3.102) 

as  in  the  previous  subsection. 

At  high  frequencies  p  =  5  =  V^  .   The  high  frequency  limiting 
value  is  found  to  be  for  any  r\ 

Re  P(co)  =  -^     .  (6.3.103) 

(jO-w       ^u 
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Fig.  7 
Noise  spectra  for  nonembedded  sphere  CH^O) 
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For  r\  =   0     equation  (6.3.91)  reduces  to  equation  (6.3.82).  The 
noise  spectrum  is  plotted  for  different  values  of  C  and  r\,   and  is 
shown  in  fig. 7.  We  distinguish  the  following  options. 

6. S. 3(a).     For  K  »   1,  or  t  »  t.,  we  have  a  surface  limited 
spectrum.  The  low  frequency  plateau  is 

Re  P(0)sa  I  . 

_2 
And,  at  high  frequencies  the  spectrum  goes  as  u   and  the  high  frequency 

limiting  value  is  found  to  be 

Re  P  (CO)  w  -^  , 
co-»«>    ^u 

as  predicted. 

6.3.3(b).      For  5  «  1  in  addition  to  r\  «  1     we  found  a 
diffusion  limited  range.  The  low  frequency  limit  is 

1  -3/2 

And  in  the  range  1  «  u  <  — =•  we  see  that  the  spectrum  goes  as  u    , 

C 
reflecting  the  typical  diffusion  behavior. 

_2 
For  high  frequencies  the  spectra  go   as  u  .  The  computer 

program  used  to  plot  the  spectrum  is  given  in  the  Appendix. 
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CHAPTER  7 
NONSYMMETRICAL  NONEMBEDDED  BODIES 


7.1  Introduction 


In  this  chapter  we  shall  analyze  the  noise  spectra  resulting 
from  diffusion  and  generation- re combination  in  nonsymmetrical  non- 
embedded  geometries.  We  have  already  seen  in  Chapter  1  the  difference 
between  embedded  and  nonembedded  bodies, the  latter  being  subjected 
to  various  boundary  conditions.   We  shall  use  the  boundary  conditions 
stemming  from  a  finite  surface  recombination  velocity  at  either  one  pair  or 
two  pairs  of  surfaces' of  the  semiconductor  sample  under' consideration. 
A  three  dimensional  case  is  of  academic  value  only, since  one  always 
has  an  ohmic  contact  in  one  of  the  directions. 

7.2   Linear  Case 
The  Green's  function  procedure  has  been  used  to  find  a  closed 
form  expression  for  the  noise  spectra  resulting  from  samples  being 
subject  to  a  finite  surface  recombination  velocity  in  the  z-direction. 

D(|^)    =   ±SAp  (7.2.1) 

±C 

where  we  have  assumed  that  the  sample  is  infinite  in  x  and  y  direction.  D  is 

the  diffusion  constant,  S  is  the  surface  recombination  velocity  and 

we  have  considered  the  sample  dimension  to  be  2C  in  the  z-direction ;Ap  being 

the  fluctuation  in  the  number  of  particles,  say  holes.  This  problem 

finds  application  in  noise  resulting  from  buried  layers  as  studied 

by  Hsieh^S.  xhe  diffusion  equation  is 
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GCz,s,z')Cs-^  7)    -   D  ^  GCz,s,z')      ^     5(z.z.).  (7.2.3) 


L(AP)     =     l^(Ap).^-^^     =     0,  (7.2.2) 

3z 

where  t  is  the  carrier  lifetime;   the  Laplace  transformed  Green's 

fimction  satisfies  the  equation 

.2_, 
GO 

8z 

According  to  the  standard  techniques  the  Green's  may  be  found 
by  setting 

,    r  uCz)  u(z')      z  5  z' 

^^^^'^'^   =  -  Dwar;^(u(z')v(z)         z  >  z'  f^-^-^) 

where  W  is  the  Wronskian  of  two  independent  solutions  u(z)  and  v(z) 
which  satisfy  the  left  and  right  side  boundary  conditions,  respectively, 
evaluated  at  z  =  z'.  The  solution  is  found  to  be  [by  Lax  and  Mengert^] 

G(z,jco,z')  =  (2YDa)   (coshyz+asinhYz) (coshyz* -asinhyz')        (7.2.5) 

for  z  5  z' 
and 

G(z,jco,z')   =   (2YDa)"  (coshyz'+asinhyz') (coshrz-asinhyz)       (7.2.6) 
for  z  >  z' 

where 

Y  =  V  C1+Jcot)/Dt  (7.2.6) 

and 


sinrC^    McoshyC 

a  = 


coshYC+  (^)  sinhyC 
Introducing  the  surface  lifetime  and  diffusion  time  by 


(7.2.7) 


Tg  =  C/S   and   t^  =  C^/D  (7.2.8) 


r.  -.-■.-    ^. 
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and  carrying  out  the  required  double  integration  the  spectrum  is 

then  found  to  be  (see  also  ^) , 

4N  P 
Sad  =  i:r-V-  Re{rf-^^  [1-  -.r ]}  C7.2.9) 

AP      N-+P-     1+JCOT  '^     r    S^ri     •       ■\      ^  ^.u    r 

0  0      •'       C — )il*Joyz)+rC   cothyC 
where 

TC  =  [-TTT^]^^^    ■  (7.2.10) 

A  critical  discussion  of  this  result  is  found  both  in  Champlin's^ 

paper  and  in  the  Lax-Mengert  paper.  The  spectrum  is  constant  at  low 

2  2 

frequencies.   If  a  range  exists  such  that  St./t  <  co  <  Tj/2t  ,  the 

-3/2 
spectrum  goes  as  oo  '  ,  whereas  at  high  frequencies  the  spectrum 

-2 
goes  as  CO  and  is  dominated  by  the  surface  lifetime  t  . 


7.5  Rectangular  Nonembedded  Bar 

This  kind  of  problem  has  been  attacked  by  Lax-Mengert^,  van 
Vliet-van  der  Ziel^^  and  Champlin^.  The  agreement  between  various 
methods  of  solutions  has  been  shown  by  van  Vliet**. 

The  equation  governing  the  flow  of  minority  carriers  (holes) 
is  the  ambipolar  transport  equation 

|_  (Ap)  =  -  ^  -  M£-V(Ap)  +  DV^(Ap)  (7.3.1) 

with  the  boundary  conditions  in  y  and  z  directions 

and  (7.3.2) 


■^iS 
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where  we  have  assumed  that  the  sample  is  infinite  in  the  x-direction  and 
that  the  noise  is  measured  between  probes  a  distance  2A  apart.  D  is 
the  diffusion  constant,  S  is  the  surface  recombination  velocity  and  we 
have  considered  the  sanple  dimension  2B  and  2C  in  y  and  z  directions 
respectively.  Further  if  no  sweep  effects  occur,  we  may  assume  that 
the  concentration  disturbances  to  be  independent  of  x.  Since  there 
is  no  field  in  y  and  z  directions  equation  (7.3.1)  becomes 

L(AP)    ,    aCAPi  ,  AP  _  DaicAPX  .  Da^   __  ^^ 

3y        dz 

We  substitute 

"  1 

L  =  S  -  -  -  A 

T 

and  C7.3.4) 

2     2 

ACy  z)    =  DC^  +  ^)   . 

9y    6z 
Note  that  A  is  slightly  differently  defined  as  in  Chapter  2.  Thus  we 
must  solve  for  the  eigenfunctions  of  A  ,  i.e.  of  the  two  dimensional 
diffusion  equation 

DC^*^)  =  Xcp   .  (7.3.5) 

3y^   az 

Subject  to  boundary  conditions  as  given  in  (7.3.2),  the 

normalized  eigenfunctions  are  found  to  be 

4V^  cos  (-g—)  cos  (-^) 


V^C2^.+sin2C^)  '"(Ti^+sin2rij) 


where  V   is  the  volume  of  the  sample  8ABC.  The  eigenvalues  are 
given  by 
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X.  .  +  DC?/B^  +  Dti?/C^  =  0  ;  (7.3.7) 

here  £,.      and  r\.     are  the  solutions  of  the  transcendental  equations 
^       J 

?tanC  =  SB/D 
and  (7.3.8) 

ritanri  =  SC/D 

So  the  Green's  function  is  given  by  (see  equation  (2.2.41)) 

cp..Cy,z)  (p..(y',z') 

G(y,z,s,y',z')   =   Z  ^^ ^ (7.3.9) 

i,j      -X. ,+S+  — 

and  the  spectrum  of  the  total  hole  fluctuation  AP  is  obtained  by  the 

standard  procedure  (by  van  Vliet  and  Fassett^) 

2     2 

'-'^  =  v^  i!,  ^\  u,i  ^^^    • 

The  normalized  spectrum  is  found  by  letting   orr  =  oo  and  multiplying 

2   -       2 
equation  (7.3.10)  by  t,  writing  co  =  Dt/B  ,  oo  =  Dt/C 

y  z 

|F.  .|^(1+  W  ?^  +  W  T)^) 

S^p(m)  =  const  2  —^ ^  J  _  ^  J  (7.3.11) 


i,j    CO  +(1+  COy^  +  W^TIj)' 


with 


B  C  (p..(y,z) 

F.  .  =   /  /  -^ dy  dz 

■'      -B  -C  a/SABC 


4      sing.  sin-Hj 


1/2,^^  _.  ^„  a/2 


Vl^  (2q+sin2q)^/''(2ri^+sin2ri^) 


C7.3.12) 
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A  particular  simple  form  is  obtained  if  S  =  «  ^  (Dirichlet 
boundary  conditions,  Ap  (surface)  =  0),  then, 

K.      =  r,   =  y(2i+l),   i=0,l,...  (7.3.13) 

and 

2...  ,^2 


1  +  D^  C2i+1)   ^  DTr^(2j-t-l)2 


AP*-  ''   "  N~+P   4  2      2     2      2 

0  0"   ^>j=°  (2i.l)(2M)[a.^.(l.  "-^^;-^^  .   P-^^J/^^  )^] 

'^      4B  4C 

(7.3.14) 
But  besides  the  remarkable  feature  that  the  low  frequency 

4 
plateau  value  is  depressed  by  a  factor  64/n  sa  0.66,  compared  to  pure 

bulk  noise,  not  much  can  be  learned  from  the  result. 

The  computer  program  for  (7.3.11)  has  been  given  in  the  Appendix, 

which  first  solves  for  the  transcendental  equation  for  C  and  i-j  and  then  / 

solves  for  the  spectrum.  The  values  of  B  and  C  were  taken  to  be  0.5  and  0.5  or 

2 
0.05  and  0.5   cm  respectively;  D  was  taken  to  be  25  cm  /sec  and  x   as 

—  X       —  ^ 
10    and  10   sec;   for  S  we  took  500  cm/sec.  The  plot  for  the  noise 

spectrum  is  shown  in  fig.  8. 

The  noise  spectrum  is  constant  at  low  frequencies  and  smoothly 

_2 
turns  over  to  f   kind  of  behavior  at  high  frequencies,  though  the 

turnover  is  slower  than  for  a  bulk  g-r  spectrum. 
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Fig.   9 
Noise  Spectrum  for  Nonembedded  Bar 
CB=0.5  cm  for  all  plots^, 
T=10'  sec. for  upper  two  plots  and  t=10~  sec. for  lower  two  plots) 


PART  C 
NOISE  DUE  TO  SURFACE  NOISE  SOURCES 
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CHAPTER  8 
NOISE  DUE  TO  STOCHASTIC  BOUNDARY  CONDITIONS 


8 . 1  Introduction 
In  this  chapter  we  investigate  the  noise  problem  with  the 
surface  Langevin  noise  source.  In  this  case  we  have  the  stochastic 
boundary  condition 

inaC^,t)  =  KOiQ,t),         r^  €  surface  (8.1.1) 

where  m  is  a  linear  surface  operator.  We  have  mentioned  earlier 
in  chapter  2  that — because  of  the  superposition  principle —  we  can 
separately  consider  the  noise  sources  associated  with  the  bulk  noise 
source  C  and  the  surface  noise  source  ^.  We  also  assume  that  these 

noise  sources  are  uncorrelated,  otherwise  there  will  be  a  contribution 

cross 
S     ,  due  to  the  correlation  of  these  noise  sources, 
a    ' 

8.2  Summary  of  Previous  Work 
During  past  years,  1/f  behavior  of  the  noise  spectrum  in  thin 
metallic  films  has  been  attributed  to  temperature  fluctuations. 
Recently  van  Vliet  et  al.^^"*  considered  ten^jerature- fluctuation  noise 
of  thin  films  supported  by  a  substrate.  They  considered  noise  due 
to  (i)  surface  heat  transfer  and  (ii)  thermal  diffusion  in  the  sample 
and  substrate.  It  was  shown  that  the  surface  source  results  in  a 
large  1/f  range,  the  level  of  which  decreases  with  increasing  heat 
conductivity  of  the  substrate.  At  high  frequencies  the  noise  goes 
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over  into  1/f  behavior.  The  thermal  diffusion  noise  is  shovm  to  result 

in  a  1/f    spectrum  at  lower  frequencies,  going  over  to  a  1/f  :   .' '; 

spectrum  at  high  frequencies.  This  noise  also  decreases  with  increasing        ■ 

heat  conductivity  of  the  substrate.  However,  for  most  practical  cases       .  • 

1/2 
the  1/f  range  is  drowned  by  the  1/f    thermal  diffusion  noise.  •:, 

In  this  chapter  we  investigate  a  similar  problem,  associated  with 
g-r  noise.  There  is  again  a  volume  and  a  surface  contribution.  The        •  ■':^,, 
aim  of  this  section  is  to  determine  (a)  whether  the  surface  contri- 
bution gives  a  1/f  range;  (b)  whether  this  noise  can  be  dominant  over 
the  bulk  g-r  noise. 


8.3   Noise  Caused  by  Surface  Generation-Recombination 

In  this  section  we  investigate  the  noise  spectrum  resulting 
from  the  noise  source  X>-     We  will  consider  three  models. 

First,  we  consider  a  bulk  semiconductor  with  surface  recom- 
bination at  both  boundaries.  The  fluctuations  in  the  whole  sanple 
are  of  interest.  Next  we  consider  a  MOSFET,  with  an  accumulation 
layer  serving  as  the  channel,  with  surface  recombination  both  on 
the  oxide  side  and  at  the  bulk  semiconductor  side.  The  noise  concerns 
here    the  carrier  population  in  the  carrier  enriched  channel  only; 
we  take  the  mobility  and  lifetime  to  be  liniform  over  the  entire  channel- 
semiconductor  region.   Finally,  we  consider  a  MOSFET  model  with  sur- 
face recombination  at  the  oxide  side  and  reflection  of  the  current  at 
the  semiconductor  side,  with  the  mobility  and  carrier  lifetime  in 
the  channel  being  different  from  those  in  the  bulk  semiconductor. 


,Kji!4^iiSl 
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8.3.1  Surface  g-r  Noise  in  a  Bulk  Semiconductor 

The  transport  problem  is  described  by  the  following  equation 

6Ap(x  t)    _   pa%(x,t)   ^    AEObll       .      gf^^t)  (8.3.1) 

3x 
with  boundary  conditions 


3x 


-   SAp  =     ^(t),  (8.3.2a) 

x=0  'x=0 


D  ^1  +   SAp  =     ^(t).  (8.3.2b) 

'^^lx=W  x=W 

For  simplicity  we  assume  the  surface  recombination  velocity  S  to  be 

equal  at  both  surfaces.   (S  is  identical  to  a,  used  in  Chapter  6). 

The  noise  due  to  the  bulk  source  ?  was  solved  for  already  in 
the  preceding  chapters;  so  we  set  it  =  0  here,  since  we  want  to 
find  the  noise  due  to  the  source  K  ,   whose  spectrum  was  given  in 
chapter  3: 

S^  =  4MaSpQ/A  (8.3.3) 

where  p„  is  the  average   concentration  of  holes,     a  =  <Ap'7/p.,  M  is 
the  modulation  factor  ?  1,   and  A  is  the  surface  area. 

To  solve   for  the  noise,   we  make  a  Fourier  transform  of  Ap  and 

of  ^ 

00 

^  \  2TTn 

Ap(x,t)      =        Z       p(co   ,x)exp(ioo  t),  oo     =  -=-  ,n=0,±l, ..  .(8. 3.4) 

n=-°° 

00 

^(t)      =        Z         £(co^)exp(iw^t).  (8.3.5) 

The   connection  between  p  and  ^  being  linear,    there  exists   a  surface 
Green's  function  H,   such  that 


m- 
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p(co^,x)      =     HCx,0,ico^)^(cOj^J.  (8.3.6) 

For  the  connection  between  the  spectra  of  Ap  and  K  we  let     T  -^  » 

and    CO    ->■  CO     =     2Trf;     then,   since     S.    (x,x' ,f)   =  lim2T<p(w  ,x)p*(co  ,x')^ 

and     S^     =     lira2T<£cco  )^*(co  )>     we  have 
K.  j_^  n  n^s' 

S^p(jc,a',f)     =     H(x,0,ico)   HCx',0,-ico)S^,  (8.3.7) 

for  any  x,x'  5  0.  The  spatially  averaged  carrier- fluctuation  noise 

for  the  san^jle  of  thickness  W  is 

.  W  W 
^An<^^^  =  ^/  /  H(x,0,ico)  H(x',0,-ico)  dx  dx'  S^       (8.3.8) 
P       W  0  0  ^ 


or 

S  " 

TT  f   H(x,0,ico)dx  ",  (8.3.9) 

0 


Thus  for  the  noise  spectrum  we  must  know  the  response  ftinction 

fCi'^)  =  tt/  H(x,0,ico)  dx.  (8.3,10) 

"  0 

To  find  the  surface  Green's  function  H,  we  define  the  Green's 

function  of  the  Laplace  transformed  transport  equation  for     s  =  ico 

(for  equation   (8.3.1  ))   as 

icoG(x,x',ia))    -  °^  GCx,x',ico)     ^     5(x-x'),  (8.3.11) 


ex2 


where 

Dico 


a 


1 

1C0+    — 
T 


subject  to  homogeneous  boundary  corresponding  to  equations  (8.3.2a) 

and  (8.3.2b) 
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D  QGCx,x',ico)    I  .  SG(x,Jc',ico)  =     0  (8.3.12a) 


ix=0 


x=0 


and 


P  aGCx^xl^io^  I  ^  sGCx,xMco)|  =     0.  (8.3.12b) 

**^  'x=W  lx=W 

The  one  dimensional  Green's   function  defined  above  satisfies  the  condition 
of  continuity  at     x  =  x'      and  the  jump  condition 


dG(x,x'  ,ico) 
9x 


^=^'*^     =     -i.  (8.3.13) 

x=x'  -  e  a  ^  ■' 


We  also  define  the  adjoint  Green's  function 

2- 
-icoG(x,x',iw)    -  a^  G(x,x',ico)     ^     5(x-x').  (8.3.14) 

9x^ 

The  regular  and  adjoint  Green's  functions  are  related  by  the  recipro- 
city relation 

G*(x,x',iu)      =     G(x',x,ico).  (8.3.15) 

The  Fourier-transform  of  equation   (8.3.1)   is 

2 
ia:^(co,x)    -   a  ^-2i^^    =     0.  (8.3.16) 

dx'^ 

According  to  Morse  and  Feshbach^^^  we  multiply  equation  (8.3.14) 
by  p  and  equation  (8.3.16)  by  G*,  then  we  subtract  these  equations 
and  integrate  the  result  over  all  x  from  0  to  W.  We  thus  obtain,  with 
the  scalar  product  being  defined  as 

W 

(f,g)   =  /  dx  fg*,  (8.3.17) 

0 


j^ 


Irp:-: 
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for    0  5  X'  5  W. 


The  Green's  theorem  for  particle  diffusion,  considering  the 
Dirichlet-Neumann  boundary  conditions  at  x=0  and  at  x=W,  is 

d^f  d^h.     ,  ^.df  ^  S.,      -rdh*   Su*i 

x=0  x=0 

Substitution  of  equation  (8.3.19)  into  equation  (8.3.18),  noticing 
the  boundary  conditions,  yields 

G'*fO  x'  ica")       1  '^ 

^  f   ^   -'    =  ^p(co,x');   0  5  X'  S  W  .  (8.3.20) 

This  gives  for  the  Fourier  coefficient  p 

p(co,x')  =  ^  G*(0,x',ico)^  .  (8.3.21) 


Further,  we  replace  x'  by  x  in  equation  (8.3.21)  and  use  the 
reciprocity  relation  given  in  equation  (8.3.15);  we  then  obtain 

p(a),x)  =  ^G(x,0,ico)^   .  (8.3.22) 

Comparing  equations    (8.3.6)    and   (8.3.22),  we  identify  the  surface 
Green's   function  in  terms   of  the  bulk  Green's  function 

H(x, 0,103)      =     ^G(x,0,iu).  (8.3.23) 


s-^ 
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Now  that  we  have  established  the  relation  between  the  surface 
and  the  bulk  Green's  function,  we  proceed  to  find  the  bulk  Green's 

function.  Consider  the  homogeneous  equation  corresponding  to  (8.3.11) 

2 

icoGCx,x',ico)  -a  ^  GCx,x',ico)  ^  ^^  (8.3.24) 

valid  for     x  ^  x' .     At     x  =  x'     we  require  continuity  of  G  and  the 
jump  condition   (8.3.13).     We  denote  the  Green's  function  for  different 
regions  as  follows 

G       for     0  <  X  <  x'  (8.3.25) 

G^       for  x'  <  X  <  W  (8.3.26) 

So  we  have 

*a  "  "•'"""  L  ■    — •'"  L 
^b 


G^  =  Asinh  f  +  Bcosh  j-  (8.3.27)  vf* 

W-x  W-x 

G,      =     Csinh  ^  +  Ecosh  ^  (8.3.28) 


where 


L     =       \E    =     /   J^,,       .  (8.3.29) 

IM  1C0+     I/t  "^  ■' 

For  G   ,  the  b.c.   at     x=0     gives 


_,,A       ,jc       B   .   ,  x^ 
D(y-coshY-  +  t-smh— ) 


or     y—     -  SB  =  0     gives 


-  S(Asinhf  +   Bcoshf)  I  =     0 

x=0  ^  '^    lx=0 


AD 
B     =     ^   >  (8.3.30) 


and  for     G,  ,    the  b.c.    at     x  =  W     gi 


ves 


£. 


K*'.'.' ' 
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-DCtcosK  ^  +  fsinh  ^)     +  SCCsinh  ~-  +  Ecosh  ^)      =  0 
^     ^    ^     ^   x=W  ^         ^   x=W 

or  ^  +  SE  =  0  or   E  =  ~   .  (8.3.31) 


Substituting  B  and  E  gives  for  G  and  G, 


G^  =  ACsinh  ^  +  |K:osh  ^) 


%  =  C(sinh-^+  _cosh-j-)    . 


(8.3.32) 
(8.3.33) 


Now,   continuity  at     x  =  x'     gives 


.   ^x'        D         ,    X' 


W-x'        D         ,    W-x" 


A(sinh-p  +  -^Y^josh  y— )    -  C(sinh  '\^     +  -^osh  ".""  )    =0        (8.3.34) 


and  the  jump  condition  (8.3.13)  gives 


X'   D 


W-x'   D 


W-x" 


A(cosh  y—  +  gj-sinh  y— )  +  C(cosh  — j +  ^y^inh  — — )     =  —.(8.3.35) 


The  system  determinant  is 


(sinh  J-  *  ^osh  — ) 


f'^^"*^  —  *  SL^i"^  — ^ 


,  .  ,  W-x'   D    ,  W-x'. 
(smh  — j —  +  ^j-cosh  — — ) 


f       ,  W-x'   D   .  .  W-x'. 
(cosh  -y—  *  ^inh  -j— ) 


,^   D   .   .  ,  W     2D    ,  W 
Ci.  _)  sxnh  y-  .  ^  cosh  -  . 


(8.3.36) 


For  A  and  C  we  find 


.   L,  .  u  W-x'   D    ,  W-x'.  ,^ 
A  =  -(sxnh  -^  +  ^osh  -j— )/D^ 

C  =  l(sinh  f  .   I^osh  f  )/D^  . 


(8.3.37) 
(8.3.38) 


*» 
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So,  finally. 


G^(x.xiico)     =     ^Csinh  —^  +  |-cosh  ^5^21)  (sinh  ^  +  ^osh  ^)      (8.3.39) 
r 

H  =  gG^Cx,0,ico)     =     J-—  Csinh  ^  +  ^  cosh  —■-)  .  (8.3.40) 


The  response  function  is 

W 
r/-  N     1  .  ,  .  ,  W-x   D    ,  W-x.  , 
f  (ico)=  ^^  f   (sinh  -J-  +  5L*^osh  -j-)    dx 

=   j^  (1  -  cosh  -  -  -g^inh  :j-)   . 

r 
The  spectrum  is,  according  to  equations  (8.3.3)  and  (8.3.10) 

Denoting  by  S  =  |f(ico)|  ,  we  find  for  the  normalized  spectrum 

2 


(8.3.41) 
(8.3.42) 


(8.3.43) 


S(co)  =  ^ 
W 


L(l  -  cosh  J-  -  -wrsinh  -r-) 
77!!  D^  ,  .  ,  W   2D   Tw" 

3  i-t 


(8.3.44) 


We  investigate  the  following  ranges. 

Ca)  At  sufficiently  low  frequencies,   |W/L|  «  1,  so  that 


..WW       u  W   , 
sinh  T-  ~  — ,    cosh  —  ~  1 


Then 


S(co) 


VX»(ico+  -) 

2 

L^j                     1 

S+2Vt)(ico+  -) 

s 

(8.3.45) 
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where,   as  in  chapter  6,  we  defined     t     =  W/S,  t,  =  W  /D,     so  that 


d     s 


(b)  At  high  frequencies   |w/l|  »  1,  so  that 


sinh 


W        ,  W     W/L 
—  '^     cosh  —  ~  e    ; 


we  then  obtain 


S(co)     ^- 

1 
W2 

(S+\/D(io5+  i)    )   \/D/Cico+  i-) 

'^^ 

S+2V'D(ico+  i)   +    c|)v4co+  i 

To                    T 

1 

(corr 

f  • 

(8.3.46) 


The  spectrum  closely  resembles  a  Lorentzian. 


8.3.2  Surface  g-r  Noise  in  a  MOSFET:  Uniform  Parameters 

The  band  diagram  for  a  MOSFET  with  inversion    layer  is  given 
in  fig.  9.  The  quantum  mechanical  carrier  density  distribution  in  a 
MOSFET  was  conputed  by  Hsing  et  al-^°  and  is  reproduced  in  fig. 10. 
This  distribution  differs  strongly  from  the  classical  distribution 
using  the  profile  in  fig. 9  and  Fermi -Dirac  statistics.  For  our  purpose 
we  shall  approximate  the  distribution  of  fig. 10  by  a  rectangle  such 
that 

n  =  n^,    0  <  X  5  Wj  (Surface  region) 
n  =  n^,    W^  <  X  5  W2.   (Bulk  region) 


'i 
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4  POTENTIAL  ENERGY 


I      Space 
h-  Charge        -• 
Region 


Fig.  9 
Energy  band  diagram  for  a  MOSFET 
(After  Hsing  et  al .  Ref.30) 
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Though  the  carrier  density  in  the  bulk  n,  is  much  less  than  the  channel 

density  n  ,  carriers  diffuse  into  the  bulk  and  undergo  bulk  recombina- 

'      s 

tion  and  surface  recombination  at  the  other  end.   In  the  present 
subsection,  we  assume  for  simplicity  that  the  surface  recombination 
velocity  in  the  oxide  and  at  the  free  surface  of  the  semiconductor 
are  the  same.  Also,  we  assume  i-i,  D,  and  t  to  be  the  same  in  the 
channel  as  in  the  bulk.  The  transport  equation  and  b.c.  are  then 
the  same  as  in  the  previous  subsection,  substituting  p  ->  n,  and  with 

Sy  =  4Man  S/A.  (8.3.47) 

The  surface  noise  may  exceed  the  bulk  noise  in  this  model  since 
n  »  n,   and  probably  M  »  1, 

For  the  surface  Green's   function  we  have  the  same  as  before; 
in  the     response  function     we  must  now  integrate  from  0  to     W^  «  W, 


W  +W     =  W.      Thus 


W 


1 


r  1  0 


(8.3.48) 


W  W 

1  r 


Using  equation  (8.3.43),  we  find  for  the  normalized  spectrum 


S(w) 


>^2   D      ^2 
L(cosh^+  g^inh  — 


cosh 


W   D 


W, 


SL-^"^  L^ 


,,   D   .  .  ,  W     2D   ,  W 
(1.  -2-^)sinh  -  .  sL^osh^ 

S  1' 


(8.3.50) 
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This   can   also  be  rewritten  as 


S(«) 


WW  W  WW  w 

2L[sinh(^  +  2r^sinh  ^  .  ^L^oshC^;^  ^   ^T^sinh  ^] 


where 


2  WW  WW 

D  1  2  1  2  2D 

(1+     -  -)  (sinh  J— cosh  j—  +   coshy— sinh  y— )    +  gY"(^) 

WW  WW 

1  2  1  2 

R  =    (cosh  -J —  cosh  -j —  +  sinh  -r—  sinh  -z — )  . 


(8.3.51) 


In  this  result  we  assume  W  »  W  .   We  discuss  the  following 


cases. 


(a)  At  sufficiently  low  frequencies 


W, 


2L 


«   1, 


W, 


« 


so  that 


WW  W 

1  1  2 

sinh  2r    ~    TT      ^""^         ^°^^  T —     ~     1 


Thus  for  the  normalized  spectrum,  noting  W  »  W^ , 


S(m) 


oy^very  low     W 


1 


WW  W 

2      1        D        1 

'^^'■L     2L       SL  2L-' 


2  W         W 

,,      D^      ,,1  2        2D, 


(8.3.52) 


(W^-  I) 


f^      D        .  ,,.,        2D. 
(U  -^-^)  (W^^  g-) 


S"L 


(8.3.53) 


constant 


.,      D     1.2         2   D 


(8.3.54) 
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2     2 
We  have  D/S  =  t  /t,;  this  gives 


SC«) 


constant 
d 


(8.3.55) 


thus 


(b)  At  intermediate  frequencies 


W, 


2L 


«  1   and 


W, 


W^    W^        W^        W^    W2/L 
sinh  -jr^    jr,        sinh  t— ~  cosh  ~  fn  e 


The  normalized  spectrum  is 


»  1, 


S(co)   = 


W  /L  W    p   W  /L  ^ 


(1+ 


D 


W^/L   2D  W2/L 
2,2^"     ^ST 


S  L 


(8.3.56) 


,-1   D  -, 


C1+ 


D 


2D 
2,2  ■"  SL-' 


S  L 


2  2  2 
For  CO  »  1/t,  the  term  D  /S  L   dominates  the  denominator  and  D/SL 

dominates  the  numerator.  One  then  finds 


:--'5 


5(03) 


(8.3.57) 


a>-»-int.  T  00 
s 


The  noise  spectrum  shows  1/f  behavior  at  intermediate  frequencies, 
In  MOSFETs  it  has  been  shown  experimentally  that  1/f  behavior  is  ob- 
served to  very  high  frequencies . 


mi'-/. 


■•!frfi4 
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(c)      At  much  higher  frequencies 


2L 


sinh 


2L 


cosh  2L 


Wj/2L 


L       ~ 


»  1, 


w. 


W^/L 


»  1,     and 
Thus 


we  obtain  for  the  spectrum 


S(oo)     sa  — r 


Lti^sl^ 


(1+ 


D 


2   2^ 
S   L'^ 


(8.3.58) 


c»>-»-very  high     D 


u;  2   2    • 
W^   CO 


(8.3.59) 


The  noise  spectrum  goes   as     — :^  at  very  high  frequencies.     The   1/co  and 


CO 


1/to       ranges  meet  at  the  break  frequency  co  ,   such  that     co     = 

s  ^  D  ...  1  ,  ,  W  W^ 
J  =     — 2"  '     providing     co^  »  -  ,     where  we  used     '^g  =  5    '  "^d  ""  D~' 


^d  1 


W 


1 


Let     D  =  25  cm^/sec     and     W     =  400  A,      see  Fig. 10.      This  yields 

0)^  =   1.6x10       rad/sec     or     f^^  Ri  25x10       hz.      Clearly  we  have   always 

that     M„  »  — .      The  above  break  frequency  is   to  all   likelihood  never 
0  T 

observable.      If       t  =   lO"     sec,   then  the  spectrum  goes  as      1/f     for 

3  10 

2x10       <     f    <     25x10      . 


Hence,  the  1/f  range  can  indeed  be  many  decades! 

8.3.3  Surface  g-r  Noise  in  a  MOSFET,  Nonuniform  Parameters 

We  now  consider  a  MOSFET  model  where  the  material  parameters 
p.,D  and  t  (carrier  lifetime)  in  the  channel  are  different  from  those 
in  the  bulk  semiconductor.  We  assume  that  in  the  bulk  semiconductor 


^^'- 
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these  parameters  are  M-' ,  D'  and  t' .  We  have 

for  0  <  X  <  W  ,  the  parameters  are  M-,  D,t 
for  W-  <  X  <  W,  the  parameters  are  |j.'  ,  D*  ,  t'  , 

The  transport  equation  is 

I-  An(x,t)  -  D^^^^^-"^  .  ^^^^^^    =  ?(x,t).  (8.3.60) 

^^  ax^  '' 

For  the  boundar>'  conditions  we  assume  in  this  subsection: 

^(t)  (8.3.61a) 

(8.3.61b) 

So  S=0  at  the  free  end  of  the  bulk,  which  means  that  the  carriers 
are  reflected  at  x=W. 

This  model  is  mathematically  identical  to  the  model  considered 
by  van  Vliet  et  al.   for  the  surface  heat  transfer  noise,  provided 
we  make  the  following  substitutions  in  their  equations  (2. 17) -(2, 19) : 

a  ^  D,  g  ->■  S,   a'  H-  D'  . 

We  also  note  that  in  our  case  a  is  different  from  that  in  their  paper, 


^1^ 

x=0 

-  S^  x=0 

"■^ 

x=W 

=      0. 

VIZ.  , 


and 


Dim                            ,             D'lco  ,-   _   ,„. 

a     =     -. — T-7—  ,  a'      =     -: r7~T  (8.3.62) 


L     =  /^     ,  I'     =  V  ^       .  (8.3.63) 

ICO  ICO 


h^'-- 
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We  find  from  their  paper  for  the  response  function   (after 
the  necessary  substitutions) 

W, 


1 


f(ia3) 


W, 


,  ,   1   D'L,    W^coshC^)-! 
tanh  -. —  +  -ppf-p  tanh  -j-j- 


DL' 


W^ 
cosh  (-j — ) 


L  ,  "l   D  ,D'L  ,  ,  ^^2   ,        D  ,  ,  "ll 
I^Ctanh  J-  .   ^)pp-  tanh  ^  .  1  .  ^  tanh  ^] 


(8.3.64) 


We  now  discuss  the  different  cases, 
(a)  At  very  low  frequencies 


«  1, 


W, 


«  1  and 


we  have 


W^     W^        W^     W^       W^ 
tanh  1 —  ^  1 — ,   tanh  -p-  i^    jT'      cosh  -: — 


For  the  normalized  spectrum  we  obtain 


5(03) 


1 

W^W^D- 
1*     DLL- 

D'W^ 
■*■      SLL' 

D     ^ 
■"  SL     L 

(8.3.65) 


5(0))       *«       1 
o>-^very  low 


(b)  At  intermediate  frequencies 
Wj   Wj        W^ 


W, 


«   1, 


w. 


(8.3.66) 


»  1     so  that 


tanh  — p  Si  —p  ,     tanh  -r-p  Ri  1 ,        cosh  -j —  si  1 


We  get  for  the  spectrum 
S(w)      9i 


1 


W  D' 
1  D' 

+ 


DL' 


W 
5L'    "^  SL   L 


(8.3.67) 
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In  the  denominator  D'/SL'   dominates.  Thus,  for  w  »  1/t, 


S(«) 


t'^co 


(8.3.68) 


9 

ir  /m 


where  t^  =  W  /D'  ^  W^/D' ,  t^  =  W/S  ;^  W2/S. 

Equation  (8.3.68)  is  of  the  same  form  as  (8.3.5),  except  that  we 
need  the  bulk  parameters. 

The  spectrum  shows  1/f  behavior  at  intermediate  frequencies 


(c)  At  high  frequencies 


WW.. 
1  2  1 

tanh  f — :^  1,  tanh  -p-fa  1,   cosh 


For  the  spectrum  we  find 


»  1, 
L 


»  1,      and  we  have 


W^/L 


S(co) 


L    ,,      LD' 

w7^^'  mJ^ 


,,      D  .D'L        ,        D 


SL'DL' 


SL 


(8.3.69) 


Further  we  have  that      |l|  «  D/S,|l' |  «  D'/SThe  spectrum  becomes 

r,2 


S(«) 


2   2 
oy>very  high  W^   00 


(8.3.70) 


This  is  the  same  as  in  subsection  8.2.2,  see  equation  (8.3.59).  At 
high  frequencies  the  spectrum  goes  as  —j   .  Again,  the  high  frequency 

CO 

range  is  near  unobtainable,  since  the  break  frequency  is  of  the 
order  of  10   hz. 

We  see  that  even  though  the  material  parameters  are  not 
uniform  in  the  model  of  this  subsection  the  asymptotic  or  limiting 
behavior  of  the  noise  spectrum  is  the  same  as  in  the  model  with  uniform 
parameters . 


m 


,.j-  -..  ^  i-- « 
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Also,  the  noise  ranges  are  apparently  insensitive  to  the  boundary 
condition  (or  to  the  value  of  the  surface  recombination  velocity)  at 
X  =  W.  The  magnitude  of  the  noise  in  the  1/f  range  in  the  present 
model  is  higher  than  in  the  uniform  parameter  model,  since 
D(volume)/D (surface)  =  M.(volume)/M.(surface)  Si  S   (see  Hsing  et  al.^^). 
The  diffusion  through  the  channel,  represented  by  the  term 
(D/SL)(W./L)   in  equation  (8.3.67)  plays  no  role  since  the  channel 
is  very  thin. 

Altogether,  the  models  presented  in  subsections  8.3.2  and  8.3.3 
give  a  very  realistic  1/f  range,  extending  over  many  decades.  This  noise 
is  due  to  the  surface  source,  which  stems  from  the  recombination- 
through  "fast"  surface  states.  This  is  in  contrast  to  McWhorter's^^ 
model,  in  which  the  1/f  noise  is  attributed  to  fluctuations  through 
tunneling  into  spatially  distributed  "slow"  surface  states.  The  latter 
can  explain  1/f  noise  for  very  low  frequencies,  whereas  in  our  model 
the  1/f  range  begins  at  cd  ^   1/t,t  being  the  channel  carrier  lifetime. 
This  is  a  drawback  compared  to     McWhorter's  model^^.  However, 
since  in  MOSFETs  the  1/f  noise  often  is  observable  for  very  high  fre- 
quencies, at  least  part  of  the  1/f  noise  could  well  be  attributed 
to  the  mechanisms  considered  in  this  chapter. 


^' 


APPENDIX 
This  Appendix  lists  the  following  computer  programs  used  to 
calculate  the  noise  spectra. 

A.l  Nonsymmetrical  embedded  bar  for  diffusion  noise  source, 
(section  5.3.1). 

For  "P"  source  and  different  geometries  the  following  changes 
should  be  made: 

A.l. a.   For  "P"  source,  remove  the  term  'TEMP'  from  numerator 
of  F(L),  just  above  statement  '10  CONTINUE',  (Section 
5.2.1). 
A.l.b.   For  square  disc  use  WIBAR  =  0.1,  W2BAR  =  0.1,  W3BAR  =  100, 

(Section  5.3.2). 
A. I.e.   For  square  bar  use  WIBAR  =  .01,  W2BAR  =  10.0,  W3BAR  =  10.0, 

(Section  5.3.3). 
A.l.d.   For  cube  use  WIBAR  =  1.0,  W2BAR  =  1.0,  W3BAR  =  1.0,  (Section 
5.3.4). 
A. 2  Nonsymmetrical  nonembedded  bar,  (Section  7.3). 
A.  3  Symmetrical  embedded  sphere  for  n  =  0,  (Section  4.4.1). 
A. 4  Symmetrical  embedded  sphere  for  n  /  0,  (Section  4.4.2). 
A. 5  Symmetrical  nonembedded  sphere  for  n  =  0,  (Section  6.3.1). 
A. 6  Symmetrical  nonembedded  sphere  for  n  5^  0,  (Section  6.3.3). 
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A.l   Nonsymmetrical  Embedded  Bar  (Diffusion  Noise  Source) 

C 

C--    TYPE  DECLARATION 

C 

DIMENSION  F(101),W(101) 
C 

C--    DEFINE  CONSTANTS 
C 

WIBAR=1.  OE-2 

W2BAR=1.  0 

W3BAR=1.  0E2 

WMIN=1.  OE-3 

WMAX=1.  0E3 

NPTS=13 

DW=(WMAX/WMIN)**(1.  0/CNPTS-l)) 
C 

C--    IF  ICNTRL  =  1,  RANGE  OF  INTEGRATION  (RI)  IS  010  RANGEl  IN 
C--  NDIV  STEPS, 

C--  =  2,  RI  IS  0  TO  RANGEl  IN  NDIV  STEPS  +  RANGEl  TO 

C—  10*RANGE1  in  NDIV  STEPS, 

C--  =  3,  RI  IS  0  TO  RANGEl  IN  NDIV  STEPS  +  RANGEl  TO 

C--  10*RANGE1  IN  NDIV  STEPS  +  10*RANGE1  TO 

C--  100*RANGEI  IN  NDIV  STEPS. 

C 

ICNTRL=2 

NDIV=20 

RANGE1=1.  0 

STEP=RANGE1/NDIV 
C 

C--    COMPUTE  THE  FREQUENCY  RESPONSE 
C 

DO  40  L=1,NPTS 

WCL)=WMIN*DW**CL-1) 

F(L)=0.  0 

ZT=-STEP 

DZT=STEP 

DO  35   IZT=1,  ICNTRL 

DZT=DZT*10,  0**CIZT-1) 

NDIV1=NDIV 

DO  30  K=1,NDIV1 

ZT=ZT+DZT 

ET=-STEP 

DET=STEP 

DO  25   IET=1,  ICNTRL 

DET=DET*10.  0**(IET-1)  •' 

DO  20  J=1,NDIV1 

ET=ET+DET  • 

ZI=-STEP 

DZI=STEP 

DO  15   IZI=1, ICNTRL 

DZI=DZI*10.  0**(IZI-1) 
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DO  10   I=1,NDIV1 
ZI=ZI+DZI 

TEMP=W1BAR*ZT**2+W2BAR*ET**2+W3BAR*ZI**2 
F(L)  =F(L)+SINCSQCZT/2)*SINSQCET/2)*SINCSQCZI/2)*TEMP/ 
$         (TEMP**2+W(L)**2) 
10       CONTINUE 


15 

ZI=ZI+DZI 

20 

CONTINUE 

25 

ET=ET+DET 

30 

CONTINUE 

35 

ZT=ZT+DZT 

40 

CONTINUE 

C-- 

NORMALIZE 

C 

NORMALIZE  ABOUT  THE  CENTER  FREQUENCY 

Fl=F(NPTS/2+l) 

DO  50  L=1,NPTS 

FCL)=F(L)/F1 

F(L)=10.0*ALOG10CFCL)) 
50       CONTINUE 
C 

C--       PRINT  THE  RESULTS 
C 

WRITE(6,60)  (I,W(I),F(I),I=1,NPTS) 
60       FORMATC10X,15,2E17,6/) 

STOP 

END 


C 
C 


^v 


FUNCTION  SINCSQCX) 

IF   CX.EQ.0.0)  SINCSQ=1.0 

IF   CX.EQ.0.0)  RETURN 

TEMP=(SIN(X)/X)**2 

SINCSQ=TEMP 

RETURN 

END 


■-^-.. 


si 
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A. 2   Nonsymmetrical  Nonembedded  Bar 


C 

C--       TYPE  DECLARATION 
C 

DIMENSION  ETA(40),ZETA(40),FSQ(40,40), 
$  OMG(101),G(101) 

REAL       K1,K2,N0 

COMPLEX    FIJ 
C 

C--       DEFINE  ABSOLUTE  CONSTANTS 
C 

Pl=3. 141592654 

B=0.5 

C=0.5 

D=25. 

S=500 

T0U=1,  0E=3 

N=20 

WMIN=1.  E-2 

WMAX=1.  E+4 

NPTS=13. 
C 

C--       DEFINE  DEPENDENT  CONSTANTS 
C 

K1=S*B/D 

K2=S*C/D 

CO=1.0/TOU 

C1=D/B**2 

C2=D/C**2 

C3=T0U 

Dw=  ovmax/wmin;)  **  ci .  o/  (npts-i)  ) 

C 

C--       SOLVE  TRANSCENDENTAL  EQUATIONS  FOR  ZETACK)  §  ETA(K),  K=l,2,3,. 

C 

CALL  S0LTRN(ZETA,K1,N) 

CALL  S0LTRN(ETA,K2,N) 
C 

C--       EVALUATE  THE  MATRIX  FSQCI,J) 
C 

DO  10,  1=1, N 

Z2=2.0*ZETA(I) 

DO  10  J=1,N 

E2=2.0*ETA(J) 

FIJ=4.0*SIN(ZETA(I))*SIN(ETA(J))/ 
$    CSQRT(CMPLX(ZETA(I)*ETA(J)*(Z2+SIN(Z2))*(E2+SIN(E2)),0.0)) 

FSQ(I,J)=CABS(FIJ)**2 
10       CONTINUE 


^ 
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C 

C--       EVALUATE  THE  FREQUENCY  RESPONSE 

C 

DO  30  K=1,NPTS 

0MG(K)=WMIN*DW**(K-1) 

SUM=0.0 

DO  20  1=1, N 

DO  20  J=1,N 

TEMP=1+C1*C3*ZETA(I)**2+C2*C3*ETA(J)**2 

SUM=SUM+FSQ(I,J)*TEMP/(OMG(K)**2+TEMP**2) 
20       CONTINUE 

G(K)=SUM 
30       CONTINUE 
C 

C--       PRINT  THE  RESULTS 
C 

WRITE(6,40)(OMG(I),G(I),I=1,NPTS) 

40       FORMAT  (5X,2E20.5/) 

STOP 

END 
C 
C 

C--       SOLTRN  SOLVES  THE  TRANSCENDENTAL  EQUATION  X*TAN(X)=C  FOR 
C--        X(K),  K=l,2,3,.. .,N. 
C 

SUBROUTINE  SOLTRN (X,C,N) 

DIMENSION  X(NhCOC2),XO(2),WS(10) 

EXTERNAL  TRAN 

PI=3.14 

EPS=1.  OE-6 

X0(l)=.25 

C0(1)=C 

IT=50 

DO  10  1=1, N 

CALL  ZSYSTMCTRAN,EPS,6,l,X0,IT,WS,CO,IERR) 

X(I)=X0(1) 
10       X0(1)=X0(1)+PI 

RETURN 

END 


Si*.; 
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A. 3  Symmetrical  Embedded  Sphere  (ri=0) 


C 

C--       TYPE  DECLARATION 

C 

DIMENSION  W(101),F(101),B(101),C(101),D(101),E(101),G(101),H(101), 
$A(100) 
C 

C--       DEFINE  CONSTANTS 
C 

WMIN=l.E-3 

WMAX=1.E5 

NPTS=17 

DIV=  (WMAX/WMIN)  **  (1 . 0/  (NPTS- 1)  ) 
C 

C--       COMPUTE  THE  FREQUENCY  RESPONSE 
C 

DO  40  L=I,NPTS 

W(L)=WMIN*DW**(L-1) 

FCL)  =  0.0 

B(L)=SQRT(2*W(L)) 

C(L)-SQRT(.5*W(L)) 

A(L)=(2.*W(L)*(SIN(B(L))+C0S(B(L)))) 

D(L)=(2.*(C0S(B(L))-SIN(B(L))))+(4.*B(L)*C0S(B(L))) 

E(L)=(2.*(C0S(C(L))+SIN(C(L))))-(2.*B(L)*C0S(C(L))) 

G (L)  =  (2  .  *  (COS  (C  (L)  )  -SIN  (C  (L)  ) ) )  +  (2  .  *B  (L)  *COS  (C  (L)  )  ) 

H(L)  =  (EXP(-B(L))*A(L)+D(L)) 

H(L)  =  -2.+2.*W(L)+H(L) 

F(L)=F(L)  +  (H(L)/(B(L)**5)) 
C 

40       CONTINUE 
C 

C--       PRINT  AND  PLOT  THE  RESULTS 
C 

WRITE  (6,60)  (I,W(I),F(I),I=1,NPTS) 
60       FORMAT(10X,15,2E17.6/) 

STOP 

END 
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A. 4   Symmetrical  Embedded  Sphere  (r\^0) 


C 

C--       TYPE  DECLARATION 

C 

DIMENSION  W(lOl) ,F(101) ,B(101) ,C(101) ,D(101) ,E(101) ,G(101) ,H(101) , 
$A(100),P(100),Q(100),C1(100),G1(100),Q1(100) 
C 

C--       DEFINE  CONSTANTS 
C 

WMIN=1.E-1 

WMAX=1.E7 

ET=.01 

NPTS=17 

DW= (WMAX/WMIN) ** (1 . 0/ (NPTS- 1) ) 
C 

C--       COMPUTE  THE  FREQUENCY  RESPONSE 
C 

DO  40  L=1,NPTS 

W(L)=WMIN*DW**CL-1) 

FCL)=0.0 

W(L)=W(L)*ET 

B  (L)  =  ( {:SQRT(2  .  )  )  *  (SQRT  (ET+SQRT CET**2+W(L)  **2)  )  )  ) 

D (L)= ( (SQRT(2 . ) 3  * (SQRT (-ET+SQRT (ET**2+W (L) **2) ) ) ) 

A(L)=1./(CB(L)**2+D(L)**2)**S) 

C(L)=((10.*B(L)**2*D(L)**3)-(5.*B(L)**4*DCL))-(D(L)**5)) 

G(L)=((2.*B(L)*D(L))+  ((2*B(L)*D(L) *EXP (-B(L) )*COS(D(L) )) 
$  (4 . *D (L) *EXP (- B(L) ) *COS  (D (L) ) ) ) 
$+(EXP(-B(L))*SIN(D(L))*(D(L)**2-B(L)**2-4.*B(L)-4.))) 

C1(L)=(-(B(L)**5-10.*B(L)**3*D(L)**2+5.*B(L]*D(L)**4)) 

G1(L)=((B(L)**2-D(L)**23+((2.*D(L)*B(L)*EXP(-B(L))*SIN(D(L))) 
$+ (4 . *D(L) *EXP(-B(L) ) *SIN(D(L) )) ) 
$+(EXP(-B(L))*C0S(D(L))*(B(L)**2-D  (L)**2+4 .+4. *B(L) )) -4 .) 

H(L)=(((B(L)**2-DCL)**2))/(3.*((B(L)**2+D(L)**2]**2))) 

Q(L)  =(C(L)*G(L)) 

Q1(L)=(C1CL)*G1(L)) 

P(L)  =  (Q(L)+Q1CL)) 

E(L)  =  (A(L)*P(L)) 

F(L)  =  F(L)+E(L)+HCL) 
40       CONTINUE 
C 

C--       PRINT  THE  RESULTS 
C 

WRITE (6,60)  (I,W(I),F(I),I=1,NPTS) 
60       FORMAT(10X,15,2E17,6/) 

STOP 

END 
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A. 5   Symmetrical  Nonembedded  Sphere  (r\=0) 


C 

C--       TYPE  DECLARATION 

C 

DIMENSION  W(lOl) ,F(101) ,6(101) ,C(101) ,D(101) ,E(101) ,H(101) , 
$AC101),A1(101),D1(101) 
C 

C--       DEFINE  CONSTANTS 
C 

WMIN=l.E-3 

WMAX=1.E3 

NPTS=13 

G=0.1 

DW= (WMAX/WMIN) ** (1 . 0/ (NPTS-1) ) 
C 

C--       COMPUTE  THE  FREQUENCY  RESPONSE 
C 

DO  40  L=1,NPTS 

W(L)=WMIN*DW**(L-1) 

F(L)=0.0 

E  (L)  =  (  (-SIN  (E  (L) )  *COSH  (E  (L)  )  )  +  (E  (L)  *COS  (E  (L)  )  *COSH  (E  (L)  )  ) 
$+  (E  (L) *SIN  (E (L) ) *SINH  (E  (L) ) ) ) 

B  (L) =  ( (COS (E (L) ) *SINH (E (L) ) )  +  (E (L) *SIN (E (L) ) *SINH (E(L) ) ) 
$- (E (L) *COS (E (L) ) *COSH(E (L) ) ) ) 

C(L)=((COS(E(L))*SINH(E(L)))  +  (E(L)*SIN(E(L))*SINH(E(L))) 
$-  (E  (L) *COS (E (L) ] *COSH (E (L) ) ) - ( (COS (E (L) *SINH(E (L) ) ) /G) ) 

D(L)=((-SIN(E(L))*COSH(E(L)))+(E(L)*COS(E(L))*COSH(E(L))) 
$+(E(L)*SIN(E(L))*SINH(E(L)))+((SIN(E(L))*COSH(E(L)))/G)) 

A1(L)=3.*((A(L)*D(L))+(B(L)*C(L))) 

D1(L)=((W(L)**2)*G*((C(L)**2)+(D(L)**2))) 

F(L)  =  F(L)  +  (A1(L)/D1(L)) 
40       CONTINUE 
C 

C--       PRINT  THE  RESULTS 
C 

WRITE(6,60)  (I,W(I)),F(I),I=1,NPTS) 
60       FORMAT(10X,15,2E17,6/) 

STOP 

END 


Ms**^.^ 


^ 

^ 


127 


A. 6   Symmetrical  Nonembcdded  Sphere  (r|^0) 


C 

C--       TYPE  DECLARATION 

C 

DIMENS10NW(99) ,F(99) ,B(99) ,Cf99) ,E(99) ,A1(99) ,H(99) ,B1 (99) , 

$F1 (99) ,D1 (99) ,A(99) , P(99) ,Q(99) ,G1(99) ,P1(99) ,C1(99) ,G1(99) 

C 

C--       DEFINE  CONSTANTS 
C 

ET=0.1 

WMIN=l.E-3 

WMAX=1.E5 

NPTS=17 

G=0.1 

DW=  (IVMAX/WMIN)  **  (1 . 0/  (NPTS-1)  ) 

C 

C--       COMPUTE  THE  FREQUENCY  RESPONSE 

C 

DO  40  L=1,NPTS 

W(L)=WMIN*DW**(L-1) 

F(L)=0.0 

H (L) = (( . 707)  * (SQRT (ET+SQRT (ET**2+W(L) **2) ) ) ) 

E(L)=((.707)  *(SQRT(-ET+SQRT(ET**2+W(L)**2)))) 

A(L)=((-SIN(E(L))*COSH(H(L)))  +  (ECL)*COS(E(L))*COSH(H(L))) 

$+  (H  (L)  *S1N  (E  (L)  )  *SINH  (H  (L) ) ) ) 
BCL)=((COS(E(L))*SINH(H(L)))  +  (E(L)*SIN(E(L))*SINH(H(L))) 

$- (H (L) *COS (E (L) ) *COSH (H (L) ) ) ) 

D (L) = ( (COS (E (L) ) *SINM (H (L) ) ) + (E (L) *SIN (E (L) ) *SINH (H (L) ) ) 
$-(II(L)*COS(E(L))*COSll(ll(L)))-(fCOS(n(I,))*SINIiril(L)))/G)) 

C(L)=((-S1N(E(L))*C0SH(H(L)))+(E(L)*C0S(E(L))*C0SH(H(L))) 
$+H(L)*SIN(E(L))*SlNH(H(L)))+((SlN(E(L))*COSH(H(L)))/G)) 

P1(L)  =  ((H(L)**2)-CE(L)**2)) 

C1(L)  =  ((P1(L)*C(L))-(2.*H(L)*E(L)*D(L))) 

G1(L)=((P1(L)*D(L))  +  (2.*H(L)*E(L)*C(L))) 

Q(L) = ( ( ( (H(L) **2) *C1 (L) )+ ( (E (L) **2) *C1 (L) ) ) **2) 

Q1(L)=((((H(L)**2)*G1(L))+((E(L)**2)*G1(L)))**2) 

D1(L)=G*(Q(L)+Q1(L)) 

A1(L)=((P1(L)*C1(L))-(2.*H(L)*E(L)*G1(L))) 

B1(L)=((P1(L)*G1(L))+(2.*H(L)*E(L)*C1(L))) 

P(L)=G*((C1(L)+A1(L))+(G1(L)*B1(L))) 

Fl (L)=-3 . * ( (A(L) *A1 (L) ) + (B (L) *B1 (L) ) ) 

F(L)=F(L)+((P(L)+F1(L))/(D1(L))) 
40       CONTINUE 
C 
C--       PRINT  THE  RESULTS 


C 


WRITE (6,60)  (I,W(I),F(I),I=1,NPTS) 
60       FORMAT(10X,15,2E17,6/) 
STOP 
END 


ifc* 
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